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Abstract 



The astrophysical S-factor for the proton weak capture on 3 He is calculated 
with correlated-hyperspherical-harmonics bound and continuum wave func- 
tions corresponding to realistic Hamiltonians consisting of the Argonne 
or Argonne vi$ two-nucleon and Urbana-VIII or Urbana-IX three-nucleon in- 
teractions. The nuclear weak charge and current operators have vector and 
axial-vector components, that include one- and many-body terms. All possi- 
ble multipole transitions connecting any of the p 3 He S- and P-wave channels 
to the 4 He bound state are considered. The S-factor at a p 3 He center-of- 
mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 
10.1 x 10 -20 keV b with the AV18/UIX Hamiltonian, a factor of ~ 4.5 larger 
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than the value adopted in the standard solar model. The P-wave transitions 
are found to be important, contributing about 40 % of the calculated S-factor. 
The energy dependence is rather weak: the AV18/UIX zero-energy S'-factor 
is 9.64 x 10~ 20 keV b, only 5 % smaller than the 10 keV result quoted above. 
The model dependence is also found to be weak: the zero-energy S'-factor is 
calculated to be 10.2 x 10~ 20 keV b with the older AV14/UVIII model, only 
6 % larger than the AV18/UIX result. Our best estimate for the S-factor at 
10 keV is therefore (10.1 ± 0.6) x 10~ 20 keV b, when the theoretical uncer- 
tainty due to the model dependence is included. This value for the calculated 
S-factor is not as large as determined in fits to the Super-Kamiokande data 
in which the hep flux normalization is free. However, the precise calculation 
of the S-factor and the consequent absolute prediction for the hep neutrino 
flux will allow much greater discrimination among proposed solar neutrino 
oscillation solutions. 

21.45,+v, 27.10.+h, 95.30.Cq 
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I. INTRODUCTION AND CONCLUSIONS 



A. Motivation 



Recently, there has been a revival of interest in the reaction 3 He(p,e + z/ e ) 4 He This 
interest has been spurred by the Super-Kamiokande collaboration measurements of the en- 
ergy spectrum of electrons recoiling from scattering with solar neutrinos [0-0 • Over most of 
the spectrum, a suppression ~ 0.5 is observed relative to the Standard Solar Model (SSM) 
predictions |K|. Above 12.5 MeV, however, there is an apparent excess of events. The hep 



process, as the proton weak capture on 3 He is known, is the only source of solar neutrinos 
with energies larger than about 14 MeV-their end-point energy is about 19 MeV. This fact 
has naturally led to questions about the reliability of calculations of the hep weak capture 
cross section, upon which is based the currently accepted SSM value for the astrophysical 
S"-factor at zero energy, 2.3 x 10 -20 keV b |TTJ. In particular, Bahcall and Krastev have 



shown JXJ] that a large enhancement, by a factor in the range 25-30, of the SSM S"-factor 
value given above would essentially fit the observed excess of recoiling electrons, in any 
of three different neutrino scenarios-uniform suppression of the 8 B flux, vacuum oscillations, 
and matter-enhanced oscillations 

The theoretical description of the hep process, as well as that of the neutron and proton 
radiative captures on 2 H, 3 H, and 3 He, constitute a challenging problem from the standpoint 
of nuclear few-body theory. Its difficulty can be appreciated by comparing the measured 
values for the cross section of thermal neutron radiative capture on 1 H, 2 H, and 3 He. Their 
respective values are: 334.2±0.5 mb @, 0.508±0.015 mb @, and 0.055±0.003 mb |L5|,[16 



Thus, in going from A=2 to 4 the cross section has dropped by almost four orders of mag- 
nitude. These processes are induced by magnetic-dipole transitions between the initial two- 
cluster state in relative S-wave and the final bound state. In fact, the inhibition of the A=3 
and 4 captures has been understood for a long time [|I7|. The 3 H and 4 He wave functions, 
denoted, respectively, with ^3 and \l/4 are, to a good approximation, eigenfunctions of the 
magnetic dipole operator /j,, namely /i 2 ^3 ~ fi p ^3 and fi z \&4 ~ 0, where /i p =2.793 n.m. is 
the proton magnetic moment (note that the experimental value of the 3 H magnetic moment 
is 2.979 n.m., while 4 He has no magnetic moment). These relations would be exact, if the 
3 H and 4 He wave functions were to consist of a symmetric S-wave term only, for example, 
\I/4 = 04(S) det[p \i,p I2, n 14- Of course, tensor components in the nuclear interac- 
tions generate significant D-state admixtures, that partially spoil this eigenstate property. 
To the extent that it is approximately satisfied, though, the matrix elements |/i z |\l/i + 2) 
and {^4\n z \^i + ^) vanish due to orthogonality between the initial and final states. This 
orthogonality argument fails in the case of the deuteron, since then 

/i 2 ^2 ^ Op -/in)02(S)Xo^O > ( L1 ) 

where Xm s an d Vm t are two-nucleon spin and isospin states, respectively. The magnetic 
dipole operator can therefore connect the large S-wave component 02 (S) of the deuteron to 
a T=l x So np state (note that the orthogonality between the latter and the deuteron follows 
from the orthogonality between their respective spin-isospin states). 

This quasi-orthogonality, while again invalid in the case of the proton weak capture on 
protons, is also responsible for inhibiting the hep process. Both these reactions are induced 
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by the Gamow- Teller operator, which differs from the (leading) isovector spin part of the 
magnetic dipole operator essentially by an isospin rotation. As a result, the hep weak 
capture and nd, pd, n 3 He, and p 3 H radiative captures are extremely sensitive to: (i) small 
components in the wave functions, particularly the D-state admixtures generated by tensor 
interactions, and (ii) many-body terms in the electro-weak current operator. For example, 
two-body current contributions provide, respectively, 50 % and over 90 % of the calculated 
pd |TB[ and n 3 He [0,0 cross sections at very low energies. 

In this respect, the hep weak capture is a particularly delicate reaction, for two additional 
reasons: firstly and most importantly, the one- and two-body current contributions are 



comparable in magnitude, but of opposite sign [[TI],g(J; secondly, two-body axial currents, 
specifically those arising from excitation of A isobars which have been shown to give the 
dominant contribution, are model dependent f20|-|22|. 



This destructive interference between one- and two-body currents also occurs in the n 3 He 
("/ien") radiative capture [|TT|,|T^], with the difference that there the leading components of 
the two-body currents are model independent, and give a much larger contribution than 
that associated with the one-body current. 

The cancellation in the hep process between the one- and two-body matrix elements has 
the effect of enhancing the importance of P-wave capture channels, which would ordinarily 
be suppressed. Indeed, one of the results of the present work is that these channels give 
about 40 % of the S'-factor calculated value. That the hep process could proceed as easily 
through P- as S-wave capture was not realized-or, at least, not sufficiently appreciated |23[- 



in all earlier studies of this reaction we are aware of, with the exception of Ref. |4|, where it 
was suggested, on the basis of a very simple one-body reaction model, that the 3 Po channel 
may be important. 



B. Previous Studies of the hep Capture 

The history of hep cross section calculations has been most recently reviewed by Bahcall 



and Krastev The first estimate of the cross section ]24] was based on the calculation of 
the overlap of an s-wave proton continuum wave function and a Is neutron wave function in 
4 He. It produced a large value for the S-factor, 630 x 10~ 20 keV b, and led to the suggestion 



by Kuzmin ]2l| that between 20 % and 50 % of the neutrinos in the high-energy end of the 



flux spectrum could originate from the hep reaction. Of course, as already discussed above 
and originally pointed out by Werntz and Brennan [EH], if the 4 He and p 3 He states are 
approximated, respectively, by (Is) 4 and (ls) 3 2s c configurations (2s c is the continuum wave 
function), and antisymmetrized in space, spin, and isospin, then the capture rate vanishes 



identically. Werntz and Brennan |26j attempted to relate the matrix element of the axial 
current occurring in the hep capture to that of the electromagnetic current occurring in the 
thermal neutron radiative capture on 3 He, and provided an upper limit for the hep S'-factor, 
3.7 x 10~ 20 keV b, based on an experimental upper limit of 100 fib for the 3 He(n,7) 4 He cross 
section known at the time. 

Werntz and Brennan assumed: (i) the validity of isospin symmetry, apart from differences 
in the neutron (in hen capture) and proton (in hep capture) continuum wave functions, which 
they related to each other via \ip p {r) /ip n {r)\ ~ C (C is the usual Gamow penetration 
factor); (ii) that two-body currents dominated both the weak and radiative captures, and 
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that their matrix elements could be put in relation to each other through an isospin rotation. 
These authors refined their earlier estimate for the hep S-factor in a later publication [23[, by 
using hard-sphere phase shifts to obtain a more realistic value for the ratio of the neutron to 
proton continuum wave functions, and by including the contributions due to P-wave capture 
channels. These refinements led to an S-factor value, 8.1 x 1CT 20 keV b, considerably larger 
than they had obtained previously. They found, though, that the P-waves only contribute 
at the 10 % level. 

Subsequent studies of the hep process also attempted to relate it to the hen radiative 
capture, but recognized the importance of D-state components in the 3 He and 4 He wave 
functions-these had been ignored in Refs. [p3p6| -, and used the Chemtob-Rho prescrip- 
tion [2lJ (with some short-range modification) for the two-body terms in the electroweak 
current operator. Tegner and Bargholtz P7| and Wervelman et al. |IB found, using a shell- 
model description of the initial and final states, that two-body current contributions do not 
dominate the capture processes, in sharp contrast with the assumptions of Refs. P3"| , P(3f and 



the later conclusions of Refs. ||14|19p0| . These two groups as well as Wolfs et al. |L5| arrived, 
nevertheless, to contradictory results, due to the different values calculated for the ratio of 
weak to electromagnetic matrix elements. Tegner and Bargholtz obtained an S-factor 
value of (17 ± 8) x 10~ 20 keV b, the spread being due to the uncertain experimental value 
of the thermal neutron capture cross section before 1983. This prediction was sharpened 
by Wolfs et al. ||15|| , who measured the hen cross section precisely. They quoted an hep 
S-factor value of (15.3 ± 4.7) x 1CT 20 keV b. Wervelman et al. |16j also measured the hen 
cross section, reporting a value of (55 ± 3) ub in excellent agreement with the Wolfs et al. 
measurement of (54 ± 6) ub, but estimated an hep S'-factor in the range (57 ± 8) x 10~ 20 
keV b. These discrepancies are presumably due to the schematic wave functions used in the 
calculations. 

In an attempt to reduce the uncertainties in the predicted values for both the radiative 
and weak capture rates, fully microscopic calculations of these reactions were performed in 
the early nineties [0,0 , based on ground- and scattering-state wave functions obtained vari- 
ationally from a realistic Hamiltonian with two- and three-nucleon interactions. The main 
part of the electromagnetic current operator (denoted as "model independent") was con- 
structed consistently from the two-nucleon interaction model. The less well known ("model 
dependent") electroweak currents associated with the excitation of intermediate A isobars 
and with transition couplings, such as the electromagnetic or axial pn current, were also 
included. However, it was emphasized that their contribution was to be viewed as numeri- 
cally uncertain, as very little empirical information is available on their coupling constants 
and short-range behavior. These studies showed that both the hen and hep reactions have 
large (in the case of the radiative capture, dominant) contributions from two-body currents. 
Indeed, the values obtained with one-body only and full currents for the hep S'-factor (radia- 
tive capture cross section) were, respectively, 5.8 x 10~ 20 and 1.3 x 10~ 20 keV b (6 and 112 
lib). These results indicated that the common practice of inferring the hep S-factor from 
the measured radiative capture cross section is bound to be misleading, because of different 
initial-state interactions in the n 3 He and p 3 He channels, and because of the large contri- 
butions associated with the two-body components of the electroweak current operator, and 
their destructive interference with the one-body current contributions. Yet, the substantial 
overprediction of the hen cross section, 112 pb versus an experimental value of 55 fib, was 
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unsatisfactory. It became clear that the contributions of the "model dependent" currents, 
particularly those due to the A isobar, were unreasonably large (about 40 pb out of the total 
112 pb). It was therefore deemed necessary to include the A degrees of freedom explicitly 
in the nuclear wave functions, rather than eliminate them in favor of effective two-body 
operators acting on nucleon coordinates, as it had been done in earlier studies. This led to 
the development of the transition-correlation operator (TCO) method ||TT|— a scaled-down 
approach to a full iV+A coupled-channel treatment. The radiative capture cross section 
was now calculated to be between 75 and 80 fib |L1] (excluding the small contribution of 
the "uncertain" unj current), the spread depending on whether the nNA coupling constant 
in the transition interactions is taken either from experiment or from the quark model. In 
this approach, the hep S'-factor was calculated to be in the range between 1.4 x 10 -20 and 
3.1 x 10 -20 keV b fTI|| , the spread due to whether the axial NA coupling was determined 
by fitting the Gamow- Teller matrix element in tritium /3-decay or, again, taken from the 
quark model (uncertainties in the values of the irN A coupling had a much smaller impact). 
In fact, the SSM value for the hep S'-factor now quoted in the literature 0,0 is the average 
of these last two results. 



C. Overview of Present Calculations 



Improvements in the modeling of two- and three-nucleon interactions and the nuclear 
weak current, and the significant progress made in the last few years in the description of 
the bound and continuum four-nucleon wave functions, have prompted us to re-examine 
the hep reaction. The nuclear Hamiltonian has been taken to consist of the Argonne fig 
two-nucleon [^8| and Urbana-IX three-nucleon [29| interactions. To make contact with 
ri],B0 , however, and to have some estimate of the model dependence 

and Urbana-VIII three-nucleon 



the earlier studies 
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of the results, the older Argonne v±4 two-nucleon 
interaction models have also been used. Both these Hamiltonians, the AV18/UIX and 
AV14/UVIII, reproduce the experimental binding energies and charge radii of the trinucleons 
and 4 He in exact Green's function Monte Carlo (GFMC) calculations [0,0 • 

The correlated-hyperspherical-harmonics (CHH) method is used here to solve variation- 
ally the bound- and scattering-state four-nucleon problem [R3,0]. The binding energy of 



4 He calculated with the CHH method [|3"4"| , |3"6f is within 1-2 %, depending on the Hamiltonian 
model, of that obtained with the GFMC method. The accuracy of the CHH method to cal- 
culate scattering states has been successfully verified in the case of the trinucleon systems, 
by comparing results for a variety of Nd scattering observables obtained by a number of 
groups using different techniques |J7|. Indeed, the numerical uncertainties in the calculation 
of the trinucleon continuum have been so drastically reduced that Nd scattering observ- 
ables can now be used to directly study the sensitivity to two- and three-nucleon interaction 



models-the A y "puzzle" constitutes an excellent example of this type of studies [[J£ 



Studies along similar lines show [[39] that the CHH solutions for the four-nucleon con- 



tinuum are also highly accurate. The CHH predictions JJIJ for the n 3 H total elastic cross 



section, ot = n ( |a s | + 3 |a t | ), and coherent scattering length, 



a s /4 + 3a t /4, mea- 



sured by neutron interferometry techniques-a s and a t are the singlet and triplet scattering 
lengths-have been found to be in excellent agreement with the corresponding experimental 
values. The n 3 H cross section is known over a rather wide energy range, and its extrapola- 
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tion to zero energy is not problematic ||0 . The situation is different for the p 3 He channel, 



for which the scattering lengths have been determined from effective range extrapolations 
of data taken above 1 MeV, and are therefore somewhat uncertain, a s = (10.8 ±2.6) fm |4T[ 
and a t = (8.1 ± 0.5) fm [[|l[] or (10.2 ± 1.5) fm p7 |. Nevertheless, the CHH results are close 



to the experimental values above. For example, the AV18/UIX Hamiltonian predicts |35 
a s = 10.1 fm and a t = 9.13 fm. 

In Refs. JTT] , pU| variational Monte Carlo (VMC) wave functions had been used to describe 
both bound and scattering states. The triplet scattering length was found to be 10.1 fm 
with the AV14/UVIII Hamiltonian model, in satisfactory agreement with the experimental 
determination and the value obtained with the more accurate CHH wave functions. However, 
the present work includes all S- and P-wave channels, namely 1 So, 3 Si, 3 Po, 1 Pi, 3 Pi, and 
3 P2, while all previous works only retained the 3 Si channel, which was thought, erroneously, 
to be the dominant one. 

The nuclear weak current consists of vector and axial-vector parts, with correspond- 
ing one-, two-, and many-body components. The weak vector current is constructed from 
the isovector part of the electromagnetic current, in accordance with the conserved- vector- 
current (CVC) hypothesis. Two-body weak vector currents have "model- independent" and 
"model-dependent" components. The model-independent terms are obtained from the 
nucleon-nucleon interaction, and by construction satisfy current conservation with it. The 
leading two-body weak vector current is the "7r-like" operator, obtained from the isospin- 
dependent spin-spin and tensor nucleon-nucleon interactions. The latter also generate an 
isovector "p-like" current, while additional isovector two-body currents arise from the isospin- 
independent and isospin-dependent central and momentum-dependent interactions. These 
currents are short-ranged, and numerically far less important than the 7r-like current. With 
the exception of the p-like current, they have been neglected in the present work. The 
model-dependent currents are purely transverse, and therefore cannot be directly linked to 
the underlying two-nucleon interaction. The present calculation includes the isovector cur- 
rents associated with excitation of A isobars which, however, are found to give a rather 
small contribution in weak- vector transitions, as compared to that due to the 7r-like current. 
The 7r-like and p-like weak vector charge operators have also been retained in the present 
study. 

The leading two- and many-body terms in the axial current, in contrast to the case of 
the weak vector (or electromagnetic) current, are those due to A-isobar excitation, which 



are treated within the TCO scheme. This scheme has in fact been extended UE3] to include 



three-body connected terms which were neglected in the earlier work ||11|| . The axial charge 
operator includes the long-range pion-exchange term Kl, required by low-energy theorems 



and the partially-conserved-axial-current relation, as well as the (expected) leading short- 
range terms constructed from the central and spin-orbit components of the nucleon-nucleon 
interaction, following a prescription due to Riska and collaborators |f44 |. 

The largest model dependence is in the weak axial current. To minimize it, the poorly 
known iV — > A transition axial coupling constant has been adjusted to reproduce the exper- 
imental value of the Gamow- Teller matrix element in tritium /3-decay. While this procedure 
is inherently model dependent, its actual model dependence is in fact very weak, as has been 
shown in Ref. [05]. The analysis carried out there could be extended to the present case. 
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D. Conclusions 



We present here a discussion of the results for the astrophysical S-factor and their im- 
plications for the Super-Kamiokande (SK) solar neutrino spectrum. 



1. Results for the S -factor 

Our results for the astrophysical S-factor, defined as 

S(E) = Ea(E) exp(4 tt a/v rel ) , (1.2) 

where cr(E) is the hep cross section at center-of-mass energy E, v Te \ is the p 3 He relative 
velocity, and a is the fine structure constant, are reported in Table [I]. By inspection of 
the table, we note that: (i) the energy dependence is rather weak: the value at 10 keV is 
only about 4 % larger than that at keV; (ii) the P-wave capture states are found to be 
important, contributing about 40 % of the calculated S-factor. However, the contributions 
from D-wave channels are expected to be very small. We have verified explicitly that they are 
indeed small in 3 D! capture, (iii) The many-body axial currents associated with A excitation 
play a crucial role in the (dominant) 3 Si capture, where they reduce the S-factor by more 
than a factor of four; thus the destructive interference between the one- and many-body 
current contributions, first obtained in Ref. |20|| , is confirmed in the present study, based 
on more accurate wave functions. The (suppressed) one-body contribution comes mostly 
from transitions involving the D-state components of the 3 He and 4 He wave functions, while 
the many-body contributions are predominantly due to transitions connecting the S-state 
in 3 He to the D-state in 4 He, or viceversa. 

It is important to stress the differences between the present and all previous studies. 
Apart from ignoring, or at least underestimating, the contribution due to P-waves, the 
latter only considered the long-wavelength form of the weak multipole operators, namely, 
their q=0 limit, where q is the magnitude of the momentum transfer. In 3 P capture, for 
example, only the Co-multipole, associated with the weak axial charge, survives in this limit, 
and the corresponding S-factor is calculated to be 2.2 x 10~ 20 keV b, including two-body 
contributions. However, when the transition induced by the longitudinal component of 
the axial current (via the L -multipole, which vanishes at q—0) is also taken into account, 
the S-factor becomes 0.82 x 10 -20 keV b, because of destructive interference between the 
Co and L matrix elements (see discussion in Sec. [11 C|) . Thus use of the long- wavelength 
approximation in the calculation of the hep cross section leads to inaccurate results. 

Finally, besides the differences listed above, the present calculation also improves that 
of Ref. [llj in a number of other important respects: firstly, it uses CHH wave functions, 



corresponding to the latest generation of realistic interactions; secondly, the model for the 
nuclear weak current has been extended to include the axial charge as well as the vector 
charge and current operators. Thirdly, the one-body operators now take into account the 
1/m 2 relativistic corrections, which had previously been neglected. In 3 Si capture, for 
example, these terms increase by 25 % the dominant (but suppressed) L\ and Ei matrix 
elements calculated with the (lowest order) Gamow- Teller operator. These improvements 
in the treatment of the one-body axial current indirectly affect also the contributions of 
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the A-excitation currents, since the NA transition axial coupling constant is determined by 
reproducing the Gamow- Teller matrix element in tritium /3-decay, as discussed in Sec. [IV E 
below. 

The chief conclusion of the present work is that the hep S-f actor is predicted to be ~ 4.5 
times larger than the value adopted in the SSM. This enhancement, while very significant, is 
smaller than that first suggested in Refs. ]I]J| , and then reconsidered by the SK collaboration 
in Ref . || . A discussion of the implications of our results for the SK solar neutrino spectrum 
is given below. 

Even though our result is inherently model dependent, it is unlikely that the model 
dependence is large enough to accommodate a drastic increase in the value obtained here. 
Indeed, calculations using Hamiltonians based on the AV18 two-nucleon interaction only and 
the older AV14/UVIII two- and three-nucleon interactions predict zero energy S-factor val- 
ues of 12.1 x 10~ 20 keV b and 10.2 x 1CT 20 keV b, respectively. It should be stressed, however, 
that the AV18 model, in contrast to the AV14/UVIII, does not reproduce the experimental 
binding energies and low-energy scattering parameters of the three- and four-nucleon sys- 
tems. The AV14/UVIII prediction is only 6 % larger than the AV18/UIX zero-energy result. 
This 6 % variation should provide a fairly realistic estimate of the theoretical uncertainty 
due to the model dependence. It would be very valuable, though, to repeat the present 
study with a Hamiltonian consisting of the CD-Bonn interaction [J5| which, in contrast to 
the AV14 and AV18 models, has strongly non-local central and tensor components. We 
would expect the CD-Bonn calculation to predict an S'-factor value close to that reported 
here, provided the axial current in that calculation were again constrained to reproduce the 
known Gamow- Teller matrix element in tritium /3-decay [[45] . 

To conclude, our best estimate for the S'-factor at 10 keV cm. energy is therefore (10.1 ± 
0.6) x 10- 20 keV b. 



2. Effect on the Super- Kamiokande Solar Neutrino Spectrum 

Super-Kamiokande (SK) detects solar neutrinos by neutrino-electron scattering. The 
energy is shared between the outgoing neutrino and scattered electron, leading to a very 
weak correlation between the incoming neutrino energy and the measured electron energy. 
The electron angle relative to the solar direction is also measured, which would in principle 
allow reconstruction of the incoming neutrino energy. However, the kinematic range of 
the angle is very forward, and is comparable to the angular resolution of the detector. 
Furthermore, event-by-event reconstruction of the neutrino energy would be prevented by 
the detector background. Above its threshold of several MeV, SK is sensitive to the 8 B 
electron neutrinos. These have a total flux of 5.15 x 10 6 cm -2 s _1 in the SSM [|ll|. While 



the flux is uncertain to about 15 %, primarily due to the nuclear-physics uncertainties in 
the 7 Be(p,7) 8 B cross section, the spectral shape is more precisely known [jj7| . 

The SK results are presented as the ratio of the measured electron spectrum to that 
expected in the SSM with no neutrino oscillations. Over most of the spectrum, this ratio 
is constant at ~ 0.5. At the highest energies, however, an excess relative to 0.5 x SSM is 
seen (though it has diminished in successive data sets). The SK 825-day data, determined 
graphically from Fig. 8 of Ref. || , are shown by the points in Fig. [I] (the error bars denote the 
combined statistical and systematic error). The excess above 12.5 MeV may be interpreted 
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as neutrino-energy dependence in the neutrino oscillation probability that is not completely 
washed out in the electron spectrum. This excess has also been interpreted as possible 
evidence for a large hep flux ]1],|3]|| (though note that the data never exceeds the full SSM 
expectation from 8 B neutrinos). In the SSM, the total hep flux is very small, 2.10 x 10 3 
cm -2 s _1 . However, its endpoint energy is higher than for the 8 B neutrinos, 19 MeV instead 
of about 14 MeV, so that the hep neutrinos may be seen at the highest energies. This is 
somewhat complicated by the energy resolution of SK, which allows 8 B events beyond their 
nominal endpoint. The ratio of the hep flux to its value in the SSM (based on the hep 
S-factor prediction of Ref. |TT|) will be denoted by a, defined as 

<? 

'-'new / 1 \ 

« = o x P osc , (1.3) 

>->SSM 

where P osc is the Ziep-neutrino suppression constant. In the present work, a = (10.1 x 
10~ 20 keV b)/(2.3 x 10 -20 keV b) = 4.4, if hep neutrino oscillations are ignored. The solid 
lines in Fig. [1] indicate the effect of various values of a on the ratio of the electron spectrum 
with both 8 B and hep to that with only 8 B (the SSM). Though some differences are expected 
in the hep spectral shape due to P-wave contributions, here we simply use the standard 
hep spectrum shape j48|. In calculating this ratio, the 8 B flux in the numerator has been 



suppressed by 0.47, the best-fit constant value for the observed suppression. If the hep 
neutrinos are suppressed by ~ 0.5, then a = 2.2. Two other arbitrary values of a (10 and 
20) are shown for comparison. As for the SK data, the results are shown as a function 
of the total electron energy in 0.5 MeV bins. The last bin, shown covering 14 - 15 MeV, 
actually extends to 20 MeV. The SK energy resolution was approximated by convolution 
with a Gaussian of energy-dependent width, chosen to match the SK LINAC calibration 
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The effects of a larger hep flux should be compared to other possible distortions of the 
ratio. The data show no excess at low energies, thus limiting the size of a neutrino magnetic 



moment contribution to the scattering [50]. The 8 B neutrino energy spectrum has recently 
been remeasured by Ortiz et al. and their spectrum is significantly larger at high energies 
than that of Ref. f4~7fl . Relative to the standard spectrum, this would cause an increase in 
the ratio at high energies comparable to the a = 4.4 case. The measured electron spectrum 
is very steep, and the fraction of events above 12.5 MeV is only ~ 1% of the total above 
threshold. Thus, an error in either the energy scale or resolution could cause an apparent 
excess of events at high energy. However, these are known precisely from the SK LINAC |49| 



calibration; an error in either could explain the data only if it were at about the 3- or 4-sigma 
level [|. 

The various neutrino oscillation solutions can be distinguished by their neutrino-energy 
dependence, though the effects on the electron spectrum are small. Generally, the ratio is 
expected to be rising at high energies, much like the effect of an increased hep flux. The 
present work predicts a = 4.4 (and a = 2.2 if the hep neutrinos oscillate). From Fig. [I], this 
effect is smaller than the distortion seen in the data or found in Refs. Jl]||,[|], where the hep 
flux was fitted as a free parameter. However, the much more important point is that this is 
an absolute prediction. Fixing the value of a will significantly improve the ability of SK to 
identify the correct oscillation solution. 

In the remainder of the paper we provide details of the calculation leading to these 
conclusions. In Sec. II we derive the hep cross section in terms of reduced matrix elements of 
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the weak current multipole operators. In Sec. Ill we discuss the calculation of the bound- and 
scattering-state wave functions with the CHH method, and summarize a number of results 
obtained for the 4 He binding energy and p 3 He elastic scattering observables, comparing 
them to experimental data. In Sec. IV we review the model for the nuclear weak current 
and charge operators, while in Sec. V we provide details about the calculation of the matrix 
elements and resulting cross section. Finally, in Sec. VI we summarize and discuss our 
results. 



II. CROSS SECTION 

In this section we sketch the derivation of the cross section for the p 3 He weak capture 
process. The center-of-mass (cm.) energies of interest are of the order of 10 keV-the Gamow- 
peak energy is 10.7 keV-and it is therefore convenient to expand the p 3 He scattering state 
into partial waves, and perform a multipole decomposition of the nuclear weak charge and 
current operators. The present study includes S- and P-wave capture channels, i.e. the 
1 So, 3 Si, 3 Po, 1 Pi, 3 Pi, and 3 P2 states in the notation 2S+1 Lj with S = 0,1, and retains 
all contributing multipoles connecting these states to the J 7r =0 + 4 He ground state. The 
relevant formulas are given in the next three subsections. Note that the x Pi and 3 Pi, and 
3 Si and 3 Dx channels are coupled. For example, a pure incoming wave will produce 
both x Pi and 3 Pi outgoing waves. The degree of mixing is significant, particularly for the 
P-waves, as discussed in Sec. 1111 C. 



A. The Transition Amplitude 

The capture process 3 He(p,e + z/ e ) 4 He is induced by the weak interaction Hamiltonian p2 



Hw = ^j dxe-^+^^rW , (2-1) 

where Gy is the Fermi coupling constant (Gy—1. 14939 10~ 5 GeV~ 2 f53|), l a is the leptonic 
weak current 

^ = u v -i a {l -75K = ( I , — 1) , (2.2) 

and j CT (x) is the hadronic weak current density. The positron and (electron) neutrino mo- 
menta and spinors are denoted, respectively, by p e and p„, and v e and u v . The Bjorken 
and Drell [54] conventions are used for the metric tensor g aT and 7-matrices. However, the 
spinors are normalized as v\v e = u\u v = 1. 

The transition amplitude in the cm. frame is then given by 

(f\H w \i) = ^|r(-q; 4 He|jt(q)|p;p 3 He) , (2.3) 

where q = p e + p^, |p;p 3 He) and | — q; 4 He) represent the p 3 He scattering state with relative 
momentum p and 4 He bound state recoiling with momentum — q, respectively, and 
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r 



(q) = J dxe^j^x) = (p(q)J(q)) . (2.4) 



The dependence of the amplitude upon the spin-projections of the proton and 3 He is un- 
derstood. It is useful to perform a partial- wave expansion of the p 3 He scattering wave 
function 

= ^ E v / 2mi i (i Sl , ^^{sj^mjj^lt , (2.5) 

L5JJ Z 

with 

7r7 Lv/ ^ _ „i-i yVi ; r>-/l- ^.L's'./.l-. 

L'S' 



<t = - i i^lW*^"' > (2-6) 



where si and S3 are the proton and 3 He spin projections, L, S, and J are the relative orbital 
angular momentum, channel spin (S'=0,1), and total angular momentum (J = L + S), 
respectively, R J is the .R-matrix in channel J, and is the Coulomb phase shift, 

a L = arg[r(L + 1 + ir,)] , (2.7) 

n = ^ • (2.8) 

Urel 

Here a is the fine-structure constant and w re i is the p 3 He relative velocity, v re \ = p/fi, \i 
being the reduced mass, \i = mm 3 /(m + m 3 ) (m and m 3 are the proton and 3 He rest masses, 
respectively). Note that has been constructed to satisfy outgoing wave boundary 

conditions, and that the spin quantization axis has been chosen to lie along p, which defines 
the z-axis. Finally, the scattering wave function ^i^ Jz as well as the 4 He wave function ^4 
are obtained variationally with the correlated-hyperspherical-harmonics (CHH) method, as 
described in Sec. |TT|. 

The transition amplitude is then written as 

(f\H w \i) = %sf^ E v / 2TTTi i (^i,^ 3 |^)(5J z , m\ jj z ) 

V ^ T Q T T _ _ 



LSJJz 

'WpK^i+s V E ; A(^|e; A -jt(q)|^ 3 JJz ) 

A=0,±1 



(2.9) 



where, with the future aim of a multipole decomposition of the weak transition operators, 
the lepton vector 1 has been expanded as 

1= E l &*> (2-10) 

A=0,±1 

with l\ = e q x ■ 1 , and 

e q o = e g3 , (2.11) 

e 3 ±i = =f-j=(e ql ± ie ?2 ) • (2.12) 

The orthonormal basis e 9 i, e q2 , e q3 is defined by e q3 = q, e q2 = p x q/|p x q|, e q i = e q2 x e q3 . 
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B. The Multipole Expansion 



Standard techniques |52| can now be used to perform the multipole expansion of the 



weak charge and current matrix elements occurring in Eq. (|2.9|) . The spin quantization axis 
is along p rather than along q. Thus, we first express the states quantized along p as linear 
combinations of those quantized along q: 

ij j 2 )p = D J>jS-^ e ^) . ( 2 - 13 ) 

Ji 

where Dj, Jz are standard rotation matrices and the angles 9 and specify the 

direction q. We then make use of the transformation properties under rotations of irreducible 
tensor operators to arrive at the following expressions: 

(*4 | P f (q) | *f+3 ^) = ^(-in-y-^D^-^-e^)^^) , (2.14) 

<*4 I e^o • j f (q) | ^1 = v^F(-i) J R J - J *Di Jzfl {-^ -6, 0) L L / J (q) , (2.15) 

(*4 I e; A • jt(q) | ^ 3 JJz ) = - v^F(-i) J (-) J - Jz ^ j2 ,- A (-0, -0, 0) 

x [XMj SJ (q) + Ej SJ (q)~\ . (2.16) 

Here A = ±1, and Cj SJ , Lj SJ , Ej SJ and Mj SJ denote the reduced matrix elements of 
the Coulomb (C), longitudinal (L), transverse electric (E) and transverse magnetic (M) 
multipole operators, explicitly given by [p2| 



CuM = J ^xp(x) ji(qx)Y Uz (£) , (2.17) 
LuM = ~y dxj(x) • Vji(qx)Y UM (k) , (2.18) 
EuM = ~ q J d* j(x) • V x j ; (gx) Y£ , (2.19) 

Af H ,(<z) = y dx j(x) ■ Uqx) Y l l , (2.20) 

where Yj* are vector spherical harmonics. 

Finally, it is useful to consider the transformation properties under parity of the multipole 
operators. The weak charge/current operators have components of both scalar/polar- vector 
(V) and pseudoscalar /axial- vector (A) character, and hence 

T llz =T llz (V)+TuM , (2-21) 

where Tu z is any of the multipole operators above. Obviously, the parity of Zth-pole V- 
operators is opposite of that of Zth-pole A-operators. The parity of Coulomb, longitudinal, 
and electric Zth-pole V-operators is (— ) l , while that of magnetic Zth-pole V-operators is 
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C. The Cross Section 



The cross section for the 3 He(p,e + z/ e ) 4 He reaction at a cm. energy E is given by 

q 2 „ „\ 1 



a(E) = J 2tt5 yAm + E 



£L, — E„ — 



2m 4 y t> re i 

2 dp e dp u 



S e S v S1S3 



(2vr) 3 (2tt) 3 



^2.22) 



where Am = m + m 3 — m 4 = 19.287 MeV (m& is the 4 He rest mass), and v Te \ is the p 3 He 
relative velocity defined above. It is convenient to write: 



4 

s e s v S 1 S Z 

where the lepton tensor L UT is defined as 



(2.23) 



L t = -£>r* = -trk(l- 75 ) 



-me) r 



2E P 



7 T (l-7s) 



2E U 



= «+ « - g° ' V e ■ V U + i ^Ve.aV^ , 

with e 0123 = —1, Vg = Pe/E e and v^ = p°/E u . The nuclear tensor iV°" r is defined as 

iV CTT = wa (* sis 3 )W T *(q; Sl s 3 ) , 



where 



W=°(q; Sl 5 3 ) = ^X i5J (q; Sl s 3 )Cf J (g) 



LSJ 



w^= 3 ( q; SlS3 ) = E x o L5J (q; ^ 3 )£f J (g) 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



LSJ 



VT =±1 (q; Sl s 3 ) = _-L Ag J (q; Sl s 3 ) [±A# S J (g) + 



LSJ 



[2.28) 



The dependence upon the direction q and proton and 3 He spin projections S\ and s 3 is 
contained in the functions X^ SJ given by 



1 1 



Xf 5J (q; SlS3 ) = ^V / 2TTli L (-i) J (-) J - Jz (-s 1 ,- S 3|5J,)(5J z ,L0|JJ,) 



^2.29) 



with A = 0, ±1. Note that the Cartesian components of the lepton and nuclear tensors 
(a, t = 1,2,3) are relative to the orthonormal basis e gl , e q2 , e q3 , defined at the end of 

Sec. ira 
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The expression for the nuclear tensor can be further simplified by making use of the 
reduction formulas for the product of rotation matrices [ 55J . In fact, it can easily be shown 
that the dependence of N aT upon the angle cos 9 = p • q can be expressed in terms of 
Legendre polynomials P n (cos9) and associated Legendre functions P™(cos#) with m = 1,2. 
However, given the large number of channels included in the present study (all S- and P- 
wave capture states), the resulting equations for N aT are not particularly illuminating, and 
will not be given here. Indeed, the calculation of the cross section, Eq. fl2.22| ), is carried out 



numerically with the techniques discussed in Sec. VB 



It is useful, though, to discuss the simple case in which only the contributions involving 
transitions from the 3 Si and 3 Po capture states are considered. In the limit q = 0, one then 
finds 



a(E) 



2 Gy 

~ m 

7T V re l 



he) iLHA)r+|£r(A)r+ 



011/ 



I ^-YllO 

°0 



(A) I 



(2.30) 



where L5 X1 (A) and £^ n (A) are the longitudinal and transverse electric axial current reduced 
matrix elements (from 3 Si capture), and Cq 10 (A) is the Coulomb axial charge reduced matrix 
element (from 3 Pq capture) at q=0. Here the "Fermi function"/o(-E') is defined as 



fo(E) 



dx x V x 2 — 1 (xq — 



x) 



(2.31) 



with xq = (Am + E)/m e . The expression in Eq. ( |2.30|) can easily be related, mutatis 
mutandis, to that given in Ref. |2(| . 

Although the q=0 approximation can appear to be adequate for the hep reaction, for 
which q < 20 MeV/c and qR ~ 0.14 or less (R being the 4 He radius), the expression for 



the cross section given in Eq. ( 2.30 ) is in fact inaccurate. To elaborate this point further, 
consider the 3 Po capture. The long-wavelength forms of the C*o(g;A) and L (q;A) multi- 
poles, associated with the axial charge and longitudinal component of the axial current, are 
constant and linear in q, respectively, as can be easily inferred from Eqs. ( |2.17| )- (|2.18| ). The 
corresponding reduced matrix elements are, to leading order in q, 



no 



(?;A) 



c + 



(2.32) 
(2.33) 



where cq = Co 10 (A) in the notation of Eq. ( 2.30|) . The 3 Po capture cross section can be 
written, in this limit, as 



2 G? 



a(E; 3 P )~-^m 5 e [f (E) \c \ 2 + h(E) m 2 e \l \ 2 - 2 f 2 (E) m e M(c* l )} 



7T V rel 



(2.34) 



When the full model for the nuclear axial charge and current is considered, the constants Co 
and Iq, at zero p 3 He relative energy, are calculated to be Co = i 0.043 fm 3 / 2 and lo = i 0.197 
fm 5 / 2 (note that they are purely imaginary at E — 0). The "Fermi functions" fo(E), fi(E), 
and f2{E), that arise after integration over the phase space, at E — have the values 
/o(0) = 2.54 x 10 6 , /i(0) = 3.61 x 10 9 , and / 2 (0) = 9.59 x 10 7 . The zero energy S-factor 
obtained by including only the term Co is 2.2 xlO -20 keV b. However, when both the Cq and 
Iq terms are retained, it becomes 0.68xl0 -20 keV b. 
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In fact, this last value is still inaccurate: when not only the leading, but also the next- 
to-leading order terms are considered in the expansion of the multipoles in powers of q (see 
Sec. |V B| ) , the S"-factor for 3 P capture increases to 0.82xlCT 20 keV b, its fully converged 
value. The conclusion of this discussion is that use of the long-wavelength approximation in 
the hep reaction leads to erroneous results. 

Similar considerations also apply to the case of 3 Si capture: at values of q different from 
zero, the transition can be induced not only by the axial current via the -E'i(A) and L\(A) 
multipoles, but also by the axial charge and vector current via the Ci(A) and and Mi(V) 
multipoles. While the contribution of Mi(V) is much smaller than that of the leading -E'i(A) 
and Li(A), the contribution of Ci(A) is relatively large, and its interference with that of 



Li(A) cannot be neglected. This point is further discussed in Sec. |VIB 



As a final remark, we note that the general expression for the cross section in Eq. ( |2.22| ) 



as follows from Eqs. ( 2.23Q — ( ^.29|) contains interference terms among the reduced matrix 



elements of multipole operators connecting different capture channels. However, these in- 
terference contributions have been found to account for less than 2 % of the total S"-factor 
at zero p 3 He cm. energy. 

III. BOUND- AND SCATTERING-STATE WAVE FUNCTIONS 

The 4 He bound-state and p 3 He scattering-state wave functions are obtained variationally 
with the correlated-hyperspherical-harmonics (CHH) method from realistic Hamiltonians 
consisting of the Argonne vx$ two-nucleon PS[ and Urbana-IX three-nucleon p3] interactions 



(the AV18/UIX model), or the older Argonne Vu two-nucleon |30| and Urbana-VIII three- 
nucleon [jri]] interactions (the AV14/UVIII model). The CHH method, as implemented 
in the calculations reported in the present work, has been developed by Viviani, Kievsky, 
and Rosati in Refs. [34,3^,56.57]. Here, it will be reviewed briefly for completeness, and a 



summary of relevant results obtained for the three- and four-nucleon bound-state properties, 
and p 3 He effective-range parameters will be presented. 

A. The CHH Method 

In the CHH approach a four-nucleon wave function \1/ is expanded as 

'Ja{xa p , YAp, Z Ap ) + 1p B (*Bp, YBp, ZBp) , (3.1) 
p 

where the amplitudes ipA and ipB correspond, respectively, to the partitions 3+1 and 2+2, 
and the index p runs over the even permutations of particles ijkl. The dependence on the 
spin-isospin variables is understood. The overall antisymmetry of the wave function \1/ is 
ensured by requiring that both ip A and ipB change sign under the exchange % ^ j. 
The Jacobi variables corresponding to the partition 3+1 are defined as 

~5t-Ap — Yj , (3-2) 

y Ap = v / 4A3(r fc - Ry) , (3.3) 
z Ap = v / 372(r ; - R ijk ) , (3.4) 
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while those corresponding to the partition 2+2 are defined as 



^■Bp *j M , 

YBp = V2(Rkl — R-ij) , 

z Bp = r l — r k , 



(3.5) 
(3.6) 
(3.7) 



where R^- (R^) and Rijk denote the cm. positions of particles ij (kl) and ijk, respectively. 
In the LS'-coupling scheme, the amplitudes ipA and ip B are expanded as 



IpA^Ap, YAp, Z Ap ) = ^ Fa >P 0*( X Api VAp, Z A p) Y£, 



(3.8) 



1p B (x Bp , YBp, ZBp) — ^a,P 4>a( X Bp, l/Bp, Z Bp ) Y^p ) 



(3.9) 



where 



Y 



a,p 



[Y ha (z A p)Y ha (yAp)] e Ye 3a (*a p ) 



SiS ns aa Sk 
[titj) Taa tk 



Si 



tl 



Sot ^ J J z 



(3.10) 



Y 



B 

a,p 



[Yi la (z Bp )Y ha (yB P )} e Y i3a (x Bp 



L a 



X 



S i S j\S aa i SkSl \s ba 

[titj] Taa [tktl] Tbc 



Sot ) J Jz 



■lai ^2eti ^3ct; ^12ai 



(3.11) 



and L a ; spin 



Here a channel a is specified by: orbital angular momenta 
angular momenta S aa , Sb a , and S a ; isospins T aa and Tf, a . The total orbital and spin angular 
momenta and cluster isospins are then coupled to the assigned JJ Z and TT Z . 

The correlation factors F ap consist of the product of pair-correlation functions, that are 



obtained from solutions of two-body Schrodinger-like equations, as discussed in Ref. p4 



These correlation factors take into account the strong state-dependent correlations induced 
by the nucleon-nucleon interaction, and improve the behavior of the wave function at small 
interparticle separations, thus accelerating the convergence of the calculated quantities with 
respect to the number of required hyperspherical harmonics basis functions, defined below. 



The radial amplitudes 0„ and 0„ are further expanded as 



</>a (XAp, VAp, Z Ap ) = ^2 



n,m 



LB 



%Bp, VBp, Z Bp ) — 22 



(P) 

^4 z Ap VAp x Ap ^ x nm\ < r2pi V^p) > 

W nm(.P) ii a i 2a h a ya ( iB i \ 
^4 z Bp iJ Bp x Bp y ^nm\ i V2pi V'ip) i 



(3.12) 
(3.13) 



where the magnitudes of the Jacobi variables have been replaced by the hyperspherical 
coordinates, i.e. the hyperradius p 
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P-y X Ap + Va p + z a p — Y x Bp + Vb p + z lp 5 (3- 14) 

which is independent of the permutation p considered, and the hyperangles appropriate for 
partitions A and B. The latter are given by 

COS 3p = x Ap /p = x Bp /p , (3.15) 

cos <p2 P = VApl (p sin 3p ) , (3. 16) 

cos 02 P = Vb p I (p sin 3p ) • (3.17) 
Finally, the hyperangle functions X" m consist of the product of Jacobi polynomials 

XZntfn) = Km (sin/5) 2m P„ X2a/3Q+ ^(cos2/?) J R: iQ+ ^ /2Q+ ^(cos2 7 ) , (3.18) 

where the indices m and n run, in principle, over all non-negative integers, K^ a = Via + 
£20 + 2m + 2, and N£ m are normalization factors |34| . 



Once the expansions for the radial amplitudes 4> A and 4> B are inserted into Eqs. 
151), the wave function \l/ can schematically be written as 



7 a,A( \ a,B( \ 

^nm \rj rya.Af „ o^ i nm KrJ nra.B ( _ r\\ 

^4 Z nm KPi* 1 ) H ^4 Z nm IP' "J 



(3.19) 



where z A (p) = u(p) and z B (p) = w(p) are yet to be determined, and the factors , with 
W = A, B, include the dependence upon the hyperradius p due to the correlation functions, 
and the angles and hyperangles, denoted collectively by fl, and are given by: 

%nm \Pi ^) = ^ Fa, P Y^ p z^ p y^ p x^ p X" m (0^, (f) 3p ) . (3.20) 



The CHH method for three-nucleon systems has been most recently reviewed in Ref. |T8 



and will not be discussed here. It leads, in essence, to wave functions having the same 
structure as in Eq. (|3.19|) with suitably defined Z(p,fl). 



B. The 3 He and 4 He Wave Functions 

The Rayleigh-Ritz variational principle 

< 5 Z V\H- E\V >= (3.21) 

is used to determine the hyperradial functions z% m (p) in Eq. ( |3.19j ) and bound state energy E. 
Carrying out the variations with respect to the functions z" m leads to a set of coupled second- 
order linear differential equations in the variable p which, after discretization, is converted 
into a generalized eigenvalue problem and solved by standard numerical techniques pi 



The present status of 3 He |5l| and 4 He |3^j36|| binding energy calculations with the CHH 



method is summarized in Tables [Tl] and [TTI|. The binding energies calculated with the CHH 
method using the AV18 or AV18/UIX Hamiltonian models are within 1.5 % of corresponding 



"exact" Green's function Monte Carlo (GFMC) results |32|], and of the experimental value 
(when the three-nucleon interaction is included). The agreement between the CHH and 
GFMC results is less satisfactory when the AV14 or AV14/UVIII models are considered, 
presumably because of slower convergence of the CHH expansions for the AV14 interaction. 
This interaction has tensor components which do not vanish at the origin. 
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C. The p 3 He Continuum Wave Functions 



The p 3 He cluster wave function $ 1+ 3 Jz , having incoming orbital angular momentum L 
and channel spin S (S = 0, 1) coupled to total angular JJ Z , is expressed as 



ip^J- = + yi^JJ. ? (3 22) 

where the term ^>c vanishes in the limit of large intercluster separations, and hence de- 
scribes the system in the region where the particles are close to each other and their mutual 
interactions are strong. The term ty I ^ SJJz describes the system in the asymptotic region, 
where intercluster interactions are negligible. It is given explicitly as: 



= \ E ® Ujkl)}s' ® Y L ,(y t )} 

^ i US' 



JJz 



X <)lL'<W + Rls.l's'VP) 0U/i) > ( 3 - 23 ) 



where yj is the distance between the proton (particle i) and 3 He (particles jkl), p is the 
magnitude of the relative momentum between the two clusters, 03 is the 3 He wave function, 
and Fl and Gl are the regular and irregular Coulomb functions, respectively. The function 
g{yi) modifies the G^pyi) at small i/i by regularizing it at the origin, and g(yi) — > 1 as ?/, > 10 
fm, thus not affecting the asymptotic behavior of ^^^ Jz . Finally, the real parameters 
R-ls L'S'ip) are ^ ne -R- m atrix elements introduced in Eq. ( |2.6| ), which determine phase shifts 
and (for coupled channels) mixing angles at the energy p 2 /(2/i) (/x is p 3 He reduced mass). 
Of course, the sum over V and S' is over all values compatible with a given J and parity. 

The "core" wave function ^>c is expanded in the same CHH basis as the bound-state 
wave function, and both the matrix elements R J LS L i S >(p) and functions z® m (p) occurring in 
the expansion of ^>c are determined by making the functional 

[RUuM\ = R J ls,l>Ap) ~ ^iS JJz \H ~E 3 - ^\^i Jz ) , (3-24) 

stationary with respect to variations in the R j LS l>s' anc ^ z nm (Kohn variational principle). 
Here E 3 = —7.72 MeV is the 3 He ground-state energy. It is important to emphasize that 
the CHH scheme, in contrast to Faddeev-Yakubovsky momentum space methods, permits 
the straightforward inclusion of Coulomb distortion effects in the p 3 He channel. 

The p 3 He singlet and triplet scattering lengths predicted by the Hamiltonian models 
considered in the present work are listed in Table fH, an d ar e found in good agreement with 
available experimental values, although these are rather poorly known. The experimental 
scattering lengths have been obtained, in fact, from effective range parametrizations of data 
taken above 1 MeV, and therefore might have large systematic uncertainties. 

The most recent determination of phase-shift and mixing-angle parameters for p 3 He 
elastic scattering has been performed in Ref. [[11|] by means of an energy- dependent phase- 
shift analysis (PSA), including almost all data measured prior 1993 (for a listing of old PSAs, 
see Ref. ||41||). New measurements are currently under way at TUNL |S9]] and Madison pD 



At low energies (E < 4 MeV) the process is dominated by scattering in L=0 and 1 waves, 
with a small contribution from L=2 waves. Therefore, the important channels are: 1 Sq, 
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3 P , 3 Si- 3 D!, 1 P 1 - 3 P 1 , 3 P 2 , 1 D 2 - 3 D 2 and 3 D 3 , ignoring channels with L > 2. The general 
trend is the following: (i) the energy dependence of the S-wave phase shifts indicates that 
the L=0 channel interaction between the p and 3 He is repulsive (mostly, due to the Pauli 
principle), while that of the four P-wave phase shifts ( 3 Po, 1 Pi, 3 Pi, and 3 P 2 ) shows that in 
these channels there is a strong attraction. Indeed, this fact has led to speculations about 
the existence of four resonant states ||61|| . (ii) The D-wave phase shifts are rather tiny, even 
at E > 2 MeV. (iii) The only mixing-angle parameter playing an important role at E < 4 
MeV is e(J T = 1"), in channel ^i- 3 ?!. 

Precise measurements have been taken at a cm. energy of 1.2 MeV, and consist in differ- 
ential cross section a (9) and proton analyzing power A y {6) |6(J data (9 is the cm. scat- 
tering angle). The theoretical predictions for <j(6), obtained from the AV18 and AV18/UIX 
interactions, are compared with the corresponding experimental data in Fig. |^. Inspection 
of the figure shows that the differential cross section calculated with the AV18/UIX model 
is in excellent agreement with the data, except at backward angles. 

By comparing, in Table |TV| , the calculated phase-shift and mixing-angle parameters with 
those extracted from the PSA at E — 1.2 MeV, one observes a qualitative agreement, 
except for the 3 Pi and 3 P 2 phase shifts which are significantly underestimated in the calcu- 
lation. The mixing-angle parameter e(l~) is found to be rather large, ~ —14°, in qualitative 
agreement with that obtained from the PSA (it is worth pointing out, however, that in the 
PSA the mixing angle was constrained to vanish at E = 0, which may be unphysical). The 
experimental error for each parameter quoted in Ref. |4]J is an average uncertainty over 



the whole energy range considered, and it is therefore only indicative. It would be very 
interesting to relate these discrepancies to the Nd A y puzzle and to specific deficiencies 
in the nuclear interaction models. A detailed study of p 3 He elastic scattering is currently 
underway and will published elsewhere p3 . 



IV. THE WEAK CHARGE AND CURRENT OPERATORS 

The nuclear weak charge and current operators have scalar /polar- vector (V) and 
pseudoscalar /axial- vector (A) components 

p±(q) = p±(q;V) + p±(q;A) , (4.1) 

j±(q)=j±(q;V)+j±(q;A) , (4.2) 

where q is the momentum transfer, q = p e + p„, and the subscripts ± denote charge raising 
(+) or lowering (-) isospin indices. Each component, in turn, consists of one-, two-, and 
many-body terms that operate on the nucleon degrees of freedom: 

p(q; a) = ^ pf ] (q; a) + ^ (q; a) + . . . , (4.3) 

i i<j 

j(q;a) = ^j i (1) (q;a) + ^jg ) (q;a) + ... , (4.4) 

i i<j 

where a=V, A and the isospin indices have been suppressed to simplify the notation. 
The one-body operators p\ and have the standard expressions obtained from a non- 
relativistic reduction of the covariant single-nucleon V and A currents, and are listed below 
for convenience. The V-charge operator is written as 



20 



pr(q;V)=p i ( i 1 I i R (q;V) + pS c (q;V) 



(4.5) 



with 



pS R (q;V)=r ii± e^ 



(4.6) 



pgL(q; v) = -i (2/ !L ^ q • (<r 4 x Pi ) e iqr * 



4m 2 



The V-current operator is expressed as 



(4.7) 



Pi , e 



iqr, 



- i—n + q x crj e iq ' r< 



(4- 



where [• • • , ■ • ■] + denotes the anticommutator, p, er, and r are the nucleon's momentum, 
Pauli spin and isospin operators, respectively, and fi v is the isovector nucleon magnetic 
moment (fi v = 4.709 n.m.). Finally, the isospin raising and lowering operators are defined 
as 



r i>± = {T i>x ±ir i>y )/2. (4.9) 

The term proportional to 1/m 2 in p J - 1 RC (q; V) is the well known |64| , |65l| spin-orbit relativis- 
ts correction. The vector charge and current operators above are simply obtained from 
the corresponding isovector electromagnetic operators by the replacement Tj^/2 — > n t ±, in 
accordance with the conserved-vector-current (CVC) hypothesis. The g-dependence of the 
nucleon's vector form factors (and, in fact, also that of the axial-vector form factors below) 
has been ignored, since the weak transition under consideration here involves very small 
momentum transfers, q < 20 MeV/c. For this same reason, the Darwin-Foldy relativistic 
correction proportional to q 2 /(8m 2 ) in Pi R c(q;V) has also been neglected. The A-charge 
operator is given, to leading order, by 

(4.10) 

+ 

while the A-current operator considered in the present work includes leading and next-to- 
leading order corrections in an expansion in powers of p/m, i.e. 

j, (1) (q; A) = j$ R (q; A) +jl 1 RC (q; A) , (4.11) 

with 

j« R (q;A) = -^r i , ± o- l e i ^ , (4.12) 



P, (1) (q;A) 



9A_ 

2 m 



Pi 



• (1) , i , :j . i 

j t ,R C (q; A ) = i^^± 



9a 



P? , 



(Ti • Pi Pi , e iq ri 



1 



■ q 



P, , e iq ■« 



2 q 



(Ti • Pi , e iqr ' 



+ i q x P i e iq ' ri 



Ti,±q<Ti ■ qe q 



2mm 



(4.13) 
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The axial coupling constant g& is taken to be [6f| 1.2654±0.0042, by averaging val- 
ues obtained, respectively, from the beta asymmetry in the decay of polarized neutrons 
(1.2626±0.0033 []S^ , |58[] ) and the half-lives of the neutron and superallowed + — > + tran- 
sitions, i.e. [2/t(0+ -> + )/ft(n) - 1]=1.2681±0.0033 @. The last term in Eq. (gig) 
is the induced pseudoscalar contribution (ra M is the muon mass), for which the coupling 
constant gp is taken as gp=-6.78 As already mentioned in Sec. |, in 3 Si capture 
matrix elements of j- ^R are suppressed. Consequently, the relativistic terms included in 
JiRC wn i cn would otherwise contribute at the percent level, give in fact a 20 % contribution 

relative to that of the leading j-^R at q=0. Among these, one would naively expect the 
induced pseudoscalar term to be dominant, due to the relatively large value of gp. This is 
not the case, however, since matrix elements of the induced pseudoscalar term scale with 
gpq 2 /(2^mmJ (< 0.014 in the g-range of interest) relative to those q- j^ R (q;A). Note 

that in the limit q=0, the expressions for p^isf R (q; V) and jj- N R (q; A) reduce to the familiar 
Fermi and Gamow- Teller operators. 

In the next five subsections we describe: i) the two-body V-current and V-charge oper- 
ators, required by the CVC hypothesis; ii) the two-body A-current and A-charge operators 
due to 7T- and p-meson exchanges, and the pir mechanism; iii) the V and A current and 
charge operators associated with excitation of A-isobar resonances, treated in perturbation 
theory and within the transition-correlation-operator method. Since the expressions for 
these operators are scattered in a number of papers [|ll],^,[^,[n]] , we collect them here for 
completeness. 



A. Two-Body Weak Vector Current Operators 



The weak vector (V) current and charge operators are derived from the corresponding 
electromagnetic operators by making use of the CVC hypothesis, which for two-body terms 
implies 



(Ti,a + T jta ) , j\f (q] 7) 



1 e 



azb jy,b(q; V ) 



(4.14) 



where ji? 2 (q;7) are the isovector (charge-conserving) two-body electromagnetic currents, 



>IJ,Z 



and a, b = x,y,z are isospin Cartesian components. A similar relation holds between the 
electromagnetic charge operators and its weak vector counterparts. The charge-raising or 
lowering weak vector current (or charge) operators are then simply obtained from the linear 
combinations 



jgL(q;V) 



jgi(q;V)±ijgi(q;V) . 



(4.15) 



The two-body electromagnetic currents have "model-independent" (MI) and "model- 
dependent" (MD) components, in the classification scheme of Riska J72]. The MI terms 
are obtained from the two-nucleon interaction, and by construction satisfy current conser- 
vation with it [7D|. Studies of the electromagnetic structure of A=2-6 nuclei, such as, for 
example, the threshold electrodisintegration of the deuteron at backward angles [73], the 
magnetic form factors of the trinucleons fl42"|, the magnetic dipole transition form factors 
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in 6 Li [[F4[j, and finally the neutron and proton radiative captures on hydrogen and he- 



lium isotopes |19| , |73| , |75| -prop ert ies in which the isovector two-body currents play a large 
role and are, in fact, essential for the satisfactory description of the experimental data- 
have shown that the leading operator is the (isovector) "7r-like" current obtained from the 
isospin-dependent spin-spin and tensor interactions. The latter also generate an isovector 
"p-like" current. There are additional MI isovector currents, which arise from the central 
and momentum-dependent interactions, but these are short-ranged and have been found to 
be numerically far less important than the 7r-like current |70|,|73|. Their contributions are 



neglected in the present study. 

Use of the CVC relation leads to the 7r-like and p-like weak vector currents below: 



( 2 ) 

j- • (ki, kj, 7rV) = i (Tj x Tj)± vp S {kj)<Ti{(Ti ■ kj) - vps(h)<Tj(cri ■ ki 



kj kj 

' p _ p 



[vpsih) ~ v PS {kj)}{(Ti ■ ki)((Tj ■ kj) 



(4.16) 



$f (ki,kj; pV) = -i(r* x Tj)± v v (kj)(Ti x (tTj x kj) -v v {ki)(Tj x (a x k { ) 



Vy(ki) - Vy(kj) 
P - P 



[(^ - kj)(a l x ki) ■ {(Tj x k 3 



+((n x k^ (Tj ■ (ki x kj) + (a-j x kj) a, ■ (ki x kj)] 



ki kj 

P-P 



[vvs(ki) - v vs {kj)} 



(4.17) 



where kj and k^ are the momenta delivered to nucleons i and j with q = kj + kj, the isospin 
operators are defined as 

(Ti X Tj)± = (Ti X Tj) x ± i (Ti X Tj) y , 

and Vps(k), Vy(k), and ^vs(^) are given by 

vps(k)=v° T (k)-2v tT (k) , 
«v(A;) = v aT (k) + v tT (k) , 
v vs (k) = v T (k) , 



(4.18) 



(4.19) 
(4.20) 
(4.21) 



with 



v T (k) = 47r / r dr j (kr)v T (r) 
Jo 
4tt ro ° 



v aT (k) 
v tT (k) 



k 2 „. 
4tt roc 



^ 2 / 

ft JO 



r 2 dr [jo(kr) — l]v crT (r) 
r 2 dr ]2(kr)v tT (r) . 



(4.22) 
(4.23) 
(4.24) 



Here v T (r), v aT (r), v tT (r) are the isospin-dependent central, spin-spin, and tensor compo- 
nents of the two-nucleon interaction (either the AV14 or AV18 in the present study). The 
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factor jo{kr) — 1 in the expression for v aT (k) ensures that its volume integral vanishes. 
Configuration-space expressions are obtained from 



tif (q; a ) 



dx. e 



lqx 



dki dkj . k :r x: !k 



(2tt) 3 (2tt) 3 



(4.25) 



where a=7rV or pV. Techniques to carry out the Fourier transforms above are discussed in 
Ref. ]7D|. 

In a one-boson-exchange (OBE) model, in which the isospin-dependent central, spin- 
spin, and tensor interactions are due to tt- and p-meson exchanges, the functions Vpg(k), 
Vy(k), and vvs{k) are simply given by 



v v (k) — ► v p (k) = 
vvs(k) ->v pS = g 



ft m 



ml k 2 + ml 

p) •< i> ' 



Am 2 k 2 + m 2 



^(i+^p) 2 m 

m 2 k 
p k 2 + m 2 p 



(4.26) 
(4.27) 
(4.28) 



where m n and m p are the meson masses, f n , g p and k p are the pseudovector tcNN, vector 
and tensor pNN coupling constants, respectively, f n (k) and f p (k) denote 7riViV and pNN 
monopole form factors, i.e. 



/«(*) 



K-m 2 a 
A 2 Q + k 2 



(4.29) 



with a=ir or p. For example, in the CD-Bonn OBE model [|(J the values for the couplings 



and cutoff masses are: fH^ = 0.075, ^/47r = 0.84, k p = 6.1, = 1.7 GeV/c, and 
A p = 1.31 GeV/c. Even though the AV14 and AV18 are not OBE models, the functions 
vps(k) and, to a less extent, vy{k) and Vys(k) projected out from their v T , v aT , and v tT 
components are quite similar to those of 7r- and p-meson exchanges in Eqs. (|4.26|) — Q4.28 



(with cutoff masses of order 1 GeV/c), as shown in Refs. [70 , 75 



Among the MD (purely transverse) isovector currents, those due to excitation of A 
isobars have been found to be the most important, particularly at low momentum transfers, 
in studies of electromagnetic structure j|2] and reactions |11[ of few-nucleon systems. Their 
contribution, however, is still relatively small when compared to that of the leading 7r-like 
current. Discussion of the weak vector currents associated with A degrees of freedom is 



deferred to Sec. IV E. 



B. Two-Body Weak Vector Charge Operators 

While the main parts of the two-body electromagnetic or weak vector current are linked 
to the form of the nucleon-nucleon interaction through the continuity equation, the most 
important two-body electromagnetic or weak vector charge operators are model dependent, 
and should be viewed as relativistic corrections. Indeed, a consistent calculation of two- 
body charge effects in nuclei would require the inclusion of relativistic effects in both the 
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interaction models and nuclear wave functions. Such a program is yet to be carried out, at 
least for systems with A > 3. 

There are nevertheless rather clear indications for the relevance of two-body electro- 
magnetic charge operators from the failure of the impulse approximation in predicting the 
deuteron tensor polarization observable [7(| , and charge form factors of the three- and four- 
nucleon systems 



\7q . The model commonly used |7lJ includes the 7T-, p-, and tu-meson 
exchange charge operators with both isoscalar and isovector components, as well as the 
(isoscalar) prey and (isovector) unr'j charge transition couplings (in addition to the single- 
nucleon Darwin- Foldy and spin-orbit relativistic corrections). The n- and p- meson exchange 
charge operators are constructed from the isospin-dependent spin-spin and tensor interac- 



tions, using the same prescription adopted for the corresponding current operators ffT 



At moderate values of momentum transfer (q < 5 fm -1 ), the contribution due to the u tt- 
like" exchange charge operator has been found to be typically an order of magnitude larger 
than that of any of the remaining two-body mechanisms and one-body relativistic correc- 



tions 42 



In the present study we retain, in addition to the one-body operator of Eq. ( ^4.5| ), only the 
"7r-like"and "p-like"weak vector charge operators. In the notation of the previous subsection, 
these are given by 



pg^k^vrV) 



m 



T 



'i,± Vps(kj) cr, ■ qcTj ■ kj + r i)± v PS (ki) cr, ■ k, t <t :i ■ q 



(4.30) 



pg^k^pV) 



m 



Tj,± v v (kj)(<Ti x q) ■ (a-j x kj) 
+n,± v v {k i ){(T j x q) • (o-j x ki) 



(4.31) 



where non-local terms from retardation effects in the meson propagators or from direct cou- 



r— (2) 

plings to the exchanged mesons have been neglected J78|,[79J. In the p\ • • (ki, kj; pV) operator 



terms proportional to powers of 1/(1 + k p ), because of the large p- meson tensor coupling 
(k p ~ 6-7), have also been neglected. Indeed, these terms have been ignored also in most 
studies of nuclear charge form factors. 



C. Two-Body Weak Axial Current Operators 

In contrast to the electromagnetic case, the axial current operator is not conserved. Its 
two-body components cannot be linked to the nucleon-nucleon interaction and, in this sense, 
should be viewed as model dependent. Among the two-body axial current operators, the 
leading term is that associated with excitation of A-isobar resonances. We again defer its 
discussion to Sec. |1V h\ . In the present section we list the two-body axial current operators 
due to 7r- and p-meson exchanges (the ttA and pA currents, respectively), and the pn- 
transition mechanism (the p7rA current). Their individual contributions have been found 
numerically far less important than those from A-excitation currents in studies of weak 



transitions involving light nuclei [[^,[I|,|C)(] . These studies P0|j45|| have also found that the 
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7rA and pA current contributions interfere destructively, making their combined contribution 
almost entirely negligible. These conclusions are confirmed in the present work. 

The 7rA, pA, and p7rA current operators were first described in a systematic way by 
Chemtob and Rho pl| . Their derivation has been given in a number of articles, including 



the original reference mentioned above and the more recent review by Towner ||22| . Their 
momentum-space expressions are given by 

9a 



j^OcktfTTA) 



2m 

9A 



(Ti X Tj)± Vn{kj) CTi X kj (Tj ■ kj 



H r j± v 7T (k j ) (q + i a-i x PA cr j ■ k, + i 

m 



(4.32) 



jg^k^pA) 



2 m 



Ti x tj)± Vp(kj) q cri ■ ((Tj x kj) + i((Tj x kj) x 

-[(Ti X ((Tj X kj)] X kj 



+— T jt ± v p (kj) (a-j x k,-) x kj - i[(Ti x (crj x k,-)] x Pj + i ^± j , (4.33) 



jg^k^k^pvrA) 



fp(ki) fn(k 



_9a ^2 / T . x T \ 

m p 1 3 kf + k"j + m 



2 



X 



1 + K p ) <Tj x kj 



iP, 



(4.34) 



where P 
Pi and 



Pi. 



Pi + p^ is the sum of the initial and final momenta of nucleon i, respectively 
and the functions v T (k) and v p (k) have already been defined in Eqs. (|4.26|) - 
( |4.27| ). Configuration-space expressions are obtained by carrying out the Fourier transforms 
in Eq. ( [4.25 ). The values used for the tiNN and pNN coupling constants and cutoff masses 

6.8 fm" 1 . 
rather 
This 



are the following: ft/ An = 0.075, g 2 p /An = 0.5, n p = 6.6, K = 4.8 fm" 1 , and A f 
The p-meson coupling constants are taken from the older Bonn OBE model |81 
than from the more recent CD-Bonn interaction |4(| (g^/ATT 



16J (^/4tt = 0.81 and n p = 6.1, 
uncertainty has in fact essentially no impact on the results reported in the present work for 
two reasons. Firstly, the contribution from j^(pA), as already mentioned above, is very 
small. Secondly, the complete two-body axial current model, including the currents due to 
A-excitation discussed below, is constrained to reproduce the Gamow- Teller matrix element 
in tritium /3-decay by appropriately tuning the value of the iVA-transition axial coupling 
g* A . Hence changes in g p and k p only require a slight readjustament of the g* A value. 

Finally, note that the replacements v 7T (k) — > vps(k) and v p (k) — > vy(k) could have been 
made in the expressions for j( 2 )(7rA) and j( 2 )(pA) above, thus eliminating the need for the 
inclusion of ad hoc form factors. While this procedure would have been more satisfactory, 
since it constrains the short-range behavior of these currents in a way consistent with that 
of the two-nucleon interaction, its impact on the present calculations would still be marginal 
for the same reasons given above. 



D. Two-Body Weak Axial Charge Operators 

The model for the weak axial charge operator adopted here includes a term of pion-range 



as well as short-range terms associated with scalar- and vector-meson exchanges [§4]]. The 
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experimental evidence for the presence of these two-body axial charge mechanisms rests 
on studies of + ^ 0~ weak transitions, such as the processes 16 N(0~,120 keV)— > 16 O(0 + ) 
and 16 O(0 + )+/i-^ 16 N(0-,120 keV)+z/ M , and first-forbidden /3-decays in the lead region 
Shell-model calculations of these transitions suggest that the effective axial charge coupling 
of a bound nucleon may be enhanced by roughly a factor of two over its free nucleon value. 
There are rather strong indications that such an enhancement can be explained by two-body 
axial charge contributions [03 . 



The pion-range operator is taken as 



Pi?(ki,k*;?rA) = -i-^2 



Ti x r 



k 2 + ml 



ki + i^ j 



(4.35) 



where f w is the pion decay constant (f 7T =93 MeV), kj is the momentum transfer to nucleon 
i, and f n (k) is the monopole form factor of Eq. ( |4.29| ) with A„.=4.8 fm _1 . The structure and 
overall strength of this operator are determined by soft pion theorem and current algebra 
arguments |3^JS^] , and should therefore be viewed as "model independent" . It can also be 
derived, however, by considering nucleon-antinucleon pair contributions with pseudoscalar 
ttN coupling. 

The short-range axial charge operators can be obtained in a "model-independent" way, 
consistently with the two-nucleon interaction model. The procedure is described in Ref. ||4 



and is similar to the one used to derive the "model- independent" electromagnetic or weak 
vector currents. Here we consider the charge operators associated only with the central 
and spin-orbit components of the interaction, since these are expected to give the largest 
contributions, after the p( 2 )(7rA) operator above. This expectation is in fact confirmed in 
the present study. The momentum-space expressions are given by 

9a 



pg^k^sA) 



2m ; 



b"t,± v s (kj) + T j:± v ST (kj)] a-i-Pi + i^j 



(4.36) 



/4 2) (k 4 ,k j; vA) 



where Pj = pi + p' { , and 



9 a 

2 m 2 
■ 9a 



-l 



4 m 2 



[n,± V V (k j ) + T jt ±V VT (kj)] [(Ti ■ Pj + i ((Ti X (Tj) ■ kj 
[Ti X Tj)± V VT (kj) (Ti-ki+l^ j , 



(4.37) 



v a (k) = Airi drr 2 j Q (kr)v a (r) , 
Jo 

with a=s, st, v, and vr. The following definitions have been introduced 



v r 



3 m 
(r) + — 



2 poo 



dr r 



v [r) = -v [r) 



m 2 f°° 
— / dr' r' 
2 



hi t\ 

v r 



v b (r') 



-/\r') 
\v b \r') 



(4.38) 



(4.39) 



where v c (r), v b (r) and v bb (r) are the isospin-independent central, spin-orbit, and (L ■ S) 2 
components of the AV14 or AV18 interactions, respectively. The definitions for v ST (r) and 
v VT (r) can be obtained from those above, by replacing the isospin-independent v c (r), v b (r) 
and v bb (r) with the isospin-dependent v CT (r), v br (r) and v bbT (r). 
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E. A-Isobar Contributions 



In this section we review the treatment of the weak current and charge operators as- 
sociated with excitation of A isobars in perturbation theory and within the context of the 
transition-correlation-operator (TCO) method fi~lf . Among the two-body axial current op- 
erators, those associated with A degrees of freedom have in fact been found to be the most 
important ones [p] ,^C[|. 

In the TCO approach, the nuclear wave function is written as 



N+A 



sna + c™) 



(4.40) 



where \1/ is the purely nucleonic component, S is the symmetrizer, and the transition oper- 
ators 11™ convert NN pairs into NA and AA pairs. The latter are defined as 



U, 



TR 



AN 



1.1 



j-AA 



(4.41) 



v na = [ u ^n {r .. )a . . S . + u MI ( nj )Sy] Tl ■ T 3 , (4.42) 
U AA = [^"'{r^Si ■ S i + u^ir^Si; 1 ] T, ■ T, . (4.43) 

Here, S, and Tj are spin- and isospin-transition operators which convert nucleon i into a 
A isobar, S-j and SfJ 1 are tensor operators in which, respectively, the Pauli spin operators 
of either particle i or j, and both particles i and j are replaced by corresponding spin- 
transition operators. The vanishes in the limit of large interparticle separations, since 
no A-components can exist asymptotically. 

In the present study the * is taken from CHH solutions of the AV14/UVIII or AV18/UIX 
Hamiltonians with nucleons only interactions, while the is obtained from two-body 
bound and low-energy scattering-state solutions of the full iV-A coupled-channel problem 
with the Argonne V28Q 10 (AV28Q) interaction, containing explicit N and A degrees of 
freedom. This aspect of the present calculations, including the validity of the approximation 
inherent to Eq. (|4.40|) , were discussed at length in the original work [II ], and have been 



reviewed more recently in Ref. [ i^j , making a further review here unnecessary. The AV28Q 



interaction provided an excellent description of the iViV database available in the early 
eighties. No attempt has been made to refit this model to the more recent and much more 
extensive Nijmegen database [|5[ . 

In the TCO scheme, the perturbation theory description of A-admixtures is equivalent 
to the replacements: 



^A,PT = ^ ^ (444) 



vMN -> NA) 



m — ira 



J 2{m — rriA) 

where the kinetic energy contributions in the denominators of Eqs. ( ^4.44| ) and ( |4.45| ) have 
been neglected (static A approximation). The transition interactions Vij(NN — > iVA) and 
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Vij(NN -> A A) have the same operator structure as Uf- A and Ufi A of Eqs. (|4.42|) and Q4.43I) , 
but with the u UTOL (r) and u tTa (r) functions replaced by, respectively, 

«—(r) = ^L^^C(x), (4.46) 



^(r) = fl + - + \ ) — C 2 (x) . (4.47) 

47T 3 \ X X z J X 



Here a = II, III, x = m n r, (ff) a = f w f*, f*f* for a = II, III, respectively, f* being the tcNA 
coupling constant, and the cutoff function C{x) = 1 — e~ Xx . In the AV28Q interaction 
f* = (6y/2/5)f n , as obtained in the quark-model, and A = 4.09. When compared to U? R , 



the perturbation theory U { - ' corresponding to Eqs. ( |4.44| ) and (|4.45|) produces NA and 



AA admixtures that are too large at short distances, and therefore leads to a substantial 



overprediction of the effects associated with A isobars in electroweak observables [11 . 

We now turn our attention to the discussion of NA and AA weak transition operators. 
The axial current and charge operators associated with excitation of A isobars are modeled 

as 

j« (q; N -> A, A) = -g\ T it± S, , (4.48) 
jf } (q; A - A, A) = -g A Q it± e iqr * , (4.49) 

and 

Pi (1) (q; N - A, A) = — T ii± ^ • Pi e iqr * (4.50) 

rriA 

Pj (1) (q;A^A,A) = --^0, ± S, Jp^e'H , (4.51) 

2 m A L J + 

where m& is the A-isobar mass, S (0) is the Pauli operator for the A spin 3/2 (isospin 3/2), 
and T^± and Q it ± are defined in analogy to Eq. ( |4.9|) . The expression for (q; A — > N, A) 
(p 4 (1) (q; A -> iV, A)) is obtained from that for j 4 (1) (q; N -> A, A) (pp (q; A^ -> A, A)) by 
replacing S, and Tj by their hermitian conjugates. The coupling constants g* A and ^ are 
not well known. In the quark-model, they are related to the axial coupling constant of 
the nucleon by the relations g* A = (6^/2/5)g A and g A = (l/5)g A . These values have often 
been used in the literature in the calculation of A-induced axial current contributions to 
weak transitions. However, given the uncertainties inherent to quark-model predictions, a 
more reliable estimate for g* A is obtained by determining its value phenomenologically in the 
following way. It is well established by now |45| that the one-body axial current of Eq. Q4.12 ) 



leads to a ~ 4 % underprediction of the measured Gamow- Teller matrix element in tritium 
/3-decay, see Table [Vj Since the contributions of A — » A axial currents (as well as those due 
to the two-body operators of Sec. |I V C|) are found to be numerically very small, as can be seen 
again from Table [V|, this 4 % discrepancy can then be used to determine g* A . Obviously, 



this procedure produces different values for g* A depending on how the A-isobar degrees of 
freedom are treated. These values are listed in Table for comparison. The g* A value 
that is determined in the context of a TCO calculation based on the AV28Q interaction, 
is about 40 % larger than the naive quark-model estimate. However, when perturbation 
theory is used for the treatment of the A isobars, the g* A value required to reproduce the 
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Gamow- Teller matrix element of tritium /3-decay is much smaller than the TCO estimate, as 
expected. Finally, the N — > A axial current derived in perturbation theory from Eqs. ([4.44 ) 
and ( |4.48| ) is, of course, identical to the expression given in Refs. p]| , [45[| . 

The N — > A and A — > A weak vector currents are modeled, consistently with the CVC 
hypothesis, as 

jf } (q; N -> A, V) = -i ^ T ii± qxS, e iq ^ , (4.52) 

m 

jf } (q; A - A, V) = -i-A- 9 ii± q x e^ r < , (4.53) 

12 m 

where the iVA-transition magnetic moment fx* is taken equal to 3 n.m., as obtained from 
an analysis of jN data in the A-resonance region ||87|| , while the value used for the A 



magnetic moment fx is 4.35 n.m. by averaging results of a soft-photon analysis of pion- 



proton bremsstrahlung data near the A ++ resonance |88 |. The contributions due to the 
weak vector currents above have been in fact found to be very small in the p 3 He capture 
process. Finally, A to A weak vector charge operators are ignored in the present study, 
since their associated contributions are expected to be negligible. 



V. CALCULATION 

The calculation of the p 3 He weak capture cross section proceeds in two steps: firstly, 
the Monte Carlo evaluation of the weak charge and current operator matrix elements, and 
the subsequent decomposition of these in terms of reduced matrix elements; secondly, the 
evaluation of the cross section by carrying out the integrations in Eq. Q2.22 ). 



A. Monte Carlo Calculation of Matrix Elements 

In a frame where the direction of the momentum transfer q also defines the quantization 
axis of the nuclear spins, the matrix element of, as an example, the weak axial (or vector) 
current has the multipole expansion 

^4|e; A -jt( q )|^ 3 ^= A ) = v ^i^ [XM^ J ( q ) + E^ SJ ( q )} , (5.1) 

with A = ±1. The expansion above is easily obtained from that in Eq. (|2.16|) , in which the 
quantization axis for the nuclear spins was taken along the direction of the relative momen- 
tum p, by setting 9=(p=0 and using Dj, Jz (0, 0, 0) = <5/^j z . Then, again as an example, the 
reduced matrix element of the axial electric dipole operator involving a transition from the 
p 3 He 3 Si state is simply given by 

^ 011 (g;A) = --L(v& 4 |e; A -jt(q;A)|< 1 3 Jz=A ) . (5.2) 

V 27T 

The problem is now reduced to the evaluation of matrix elements of the same type as on 
the right-hand-side of Eq. (|5.2j ). These can schematically be written as 
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(^4,iV+A | O | ^l+3,Af+A) /g ^ 

[(^4,7V+A |^4,iV+A)(^l+3,iV+A | ^1+3,AT+a)] ^ 

where the initial and final states have the form of Eq. ( |4.40| ). It is convenient to expand the 
latter as 

^V+A = ^ + XX R * + --- ' ( 5 - 4 ) 

i<j 

so that the numerator of Eq. (|5.3|) can be expressed as 

(^4,v+a | O | ^i+ 3 ,iv+A) = (*4 I O(iVonly) | * 1+3 ) + (* 4 | 0(A) | * 1+3 ) , (5.5) 

where the operator 0(iV only) denotes all one- and two-body contributions to the weak 
charge or current operator O, involving only nucleon degrees of freedom, i.e. 0(iV only) = 
(1) (iV -> N) + 0®(NN -> NN), while 0(A) includes terms that involve the A-isobar 
degrees of freedom, associated with the explicit A transitions 0^'(N — ► A), 0^(A — > N), 
0^(A — > A), and with the transition operators U? R . A diagrammatical illustration of the 
terms contributing to 0(A) is given in Fig. ||: the terms (a)-(e), (f)-(i), and (j) represent, 
respectively, two-, three-, and four-body operators. The terms (e) and (g)-(j) are to be 
viewed as renormalization corrections to the "nucleonic" matrix element of 0^(N — > N), 
due to the presence of A-admixtures in the wave functions. Connected three-body terms 
containing more than a single A isobar have been ignored, since their contributions are 
expected to be negligible. Indeed, the contribution from diagram (d) has already been 
found numerically very small. 

The two-body terms of Fig. |3| are expanded as operators acting on the nucleons' coordi- 
nates. For example, the terms (a) and (c) in Fig. |^ have the structure, respectively, 

(a) = ^ t 0f ) (A^A), (5.6) 

{c) = u ^ w {a ^ a)u an ^ (5 7) 

which can be reduced to operators involving only Pauli spin and isospin matrices by using 
the identities: 

S t -AS-B = -A-B--o--(AxB) , (5.8) 
3 3 

S t -AE-BS-C = -iA-(BxC)--<r-AB-C 

3 v ; 3 

— ttA • B C ■ cr + ^A • (B ■ <t)C , (5.9) 

where A, B and C are vector operators that commute with cr, but not necessarily among 
themselves. While the three- and four-body terms in Fig. |3| could have been reduced in 
precisely the same way, the resulting expressions in terms of cr and r matrices become too 
cumbersome. Thus, for these it was found to be more convenient to retain the explicit 
representation of S (S" 1 ") ClS db 4 x 2 (2 x 4) matrix 
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and of S as a 4 x 4 matrix 



3e 
-\/6e H 





\/6e_ 
e V^e- 
-y/8e + — e v^e- 
—\/6e + — 3e 

(— ) M e_ M , and derive the result of terms 
A) Ujp, N on state \"$?) by first operating with U^. N , then with 

A), and finally with Uj^^ . These terms (as well as three-body contributions to 
the wave function normalizations, see below) were neglected in the calculations reported in 
Ref. 0. 

Of course, the presence of A-admixtures also influences the normalization of the wave 
functions, as is obvious from Eq. (|5.3|): 



where e± = =f(x ± iy)/\/2, e = z, and e* 
such as (f)=^ At O] 1) (A 
Of (A 



N+A | ^ N+A 



i<3 . fc 7^ ,3 



i<3 

[u% N] u? k N + u$*u»*] I *> + ... 



(5.10) 



The wave function normalization ratios (^n+a \ ^n+a)/{^ \ obtained for the bound 



three- and four-nucleon systems, are listed in Table [VI 4 Note that the normalization of the 
p 3 He continuum state is the same as that of 3 He, up to corrections of order (volume) -1 . 



The matrix elements in Eqs. (5J5) and ( 5.1C ) are computed, without any approximation, 
by Monte Carlo integrations. The wave functions are written as vectors in the spin-isospin 
space of the four nucleons for any given spatial configuration R = (ri,...,r4). For the 
given R we calculate the state vector [0(R, iVonly) + 0(R, A)] ^(R) with the techniques 
developed in Refs. p^ , [7D|] . The spatial integrations are carried out with the Monte Carlo 
method by sampling R configurations according to the Metropolis et al. algorithm [89 
using a probability density W (R) proportional to 



W(R)oc V<*I(R)*4(R)) , 



(5.11) 



where the notation (• ■ ■) implies sums over the spin-isospin states of the 4 He wave function. 
Typically 200,000 configurations are enough to achieve a relative error < 5 % on the total 
S- factor. 



B. Calculation of Cross Section 

Once the reduced matrix elements (RMEs) have been obtained, the calculation of the 
cross section cr(E) is reduced to performing the integrations over the electron and neutrino 
momenta in Eq. ( p. 22] ) numerically. We write 
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1 Cr 2 PC 



doc f. 



dx,, 



2 IT 



plf-'L aT N° 



(5.12) 



where one of the azimuthal integrations has been carried out, since the integrand only 
depends on the difference (f) = <p e — <p v . The 5-function occurring in Eq. Q2.22|) has also been 
integrated out resulting in the factor with 



./ 



Y _|_ Pe %ev Pp 



m 4 m 4 

The magnitude of the neutrino momentum is fixed by energy conservation to be 

2A 



Pv = , _ = 

l+p e x ev /m 4 + y (1 +p e x ev /m i ) 2 + 2 A/m 4 
where A = Am + E — E e — Pg/2m 4 . The variable x ev is defined as 

X ev = p e ■ p u = X e X v + y/l - xfyi - X 2 V COS , 



(5.13) 



(5.14) 



(5.15) 



where x e = cos 6 e and x v = cos 9 U . Finally, the integration over the magnitude of the electron 
momentum extends from zero up to 



Pa 



m\ + m 2 e + 2 77i4 (Am + E) 



1 2 



mi 



(5.16) 



The lepton tensor is explicitly given by Eq. ( |2.24j ), while the nuclear tensor is constructed 
using Eqs. (|2.25| )- (|2.29|) . Computer codes have been developed to calculate the required 
rotation matrices corresponding to the q-direction (9, <p) with 



cos 9 



z • q 



z ■ (p e + p„) 

|Pe + Pu\ 

Pe ■Ee Pv ■Ev 

\/pi +pI + 2p e p v x ev 



(5.17) 



Finally, note that the nuclear tensor requires the values of the RMEs at the momentum 
transfer q (the denominator in the second line of Eq. ( |5.17| ) ) . It has been found convenient 
to make the dependence upon q of the RMEs explicit by expanding 



Tj^(g) = f 



E+LSJ n 2n 
b 2n H 



(5.18) 



n>0 



consistently with Eqs. (|2.17 )- (|2.20|) . Here m = J, J± 1, depending on the RME considered. 
For example, m = 1 for the Lq 10 (A) RME. Given the low momentum transfers involved, q < 
20 MeV/c, the leading and next-to-leading order terms to and £2 are sufficient to reproduce 
accurately T(q). Note that the long- wavelength-approximation corresponds, typically, to 
retaining only the to term. 

A moderate number of Gauss points (of the order of 10) for each of the integrations in 
Eq. Q5.12D is sufficient to achieve convergence within better than one part in 10 3 . The com- 
puter program has been successfully tested by reproducing the result obtained analytically 
by retaining only the 3 Si ^i(A) and Li(A) and 3 P C (A) RMEs. 
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VI. RESULTS AND DISCUSSION 



The 5-factor calculated values are listed in Table |, and their implications to the recoil 
electron spectrum measured in the SK experiment, see Fig. [1|, have already been discussed 
in the introduction. In Tables [VII]] , FX , and KII| - [XV| , we present our results, obtained with 
the AV18/UIX Hamiltonian model, for the reduced matrix elements (RMEs) connecting any 
of the p 3 He S- and P-wave channels to the 4 He bound state. The values for these RMEs are 
given at zero energy and a lepton momentum transfer g=19.2 MeV/c. Note that the RMEs 
listed in all tables are related to those defined in Eqs. ( |2.14| )-( [2.16| ) via 



T 



-LS.J 



Ana 



[exp(47ra/f 



reW 



1] T^ SJ 



(6.1) 



which can be shown to remain finite in the limit v Te \ —>■ 0, corresponding to zero energy. 

In Table [X-V1| we list the individual contributions of the S- and P-wave capture channels 
to the total .S-factor at zero cm. energy, obtained with the AV18/UIX, the AV18 only (to 
study the effects of the three-nucleon interaction), and the older AV14/UVIII (to study the 
model dependence and to make contact with the earlier calculations of Refs. p!T|j20| ). The 
model dependence is discussed in Sec. |V1D. 



In Tables |, [VIII| , |X , and pQ| - KV| , the cumulative nucleonic contributions are normalized 



as 



[one— body +mesonic] 



(^ 4 |0(Aonly)|^ 1+3 ) 

nl/2 



(6.2) 



[(^I^X^i+al^i+s)] 

However, when the A-isobar contributions are added to the cumulative sum, the normaliza- 
tion changes to 



[one — body+ mesonic + A] 



(^4,jv+A|Q(iVonly) + Q(A)\V 1+3 , N+A ) 

1/2 



4,iV+ A I ^4, JV+ A 



l+3,iV+A | ^ l+3,iV+A 



01 



(6.3) 



As already pointed out earlier in Sec. |V Aj , the normalization of the initial scattering state is 
the same as that of 3 He, up to corrections of order (volume) -1 . In Table we also report 
results in which the A-components in the nuclear wave functions are treated in perturbation 
theory, as discussed in Sees. [IV E| and |V A| , and the 0(A) only includes the operators in 
panel (a) of Fig. |3|. In this case, the cumulative contributions [one-body+mesonic+Apx] are 
normalized as in Eq. 



A. 1 So Capture 

The 1 So capture is induced by the weak vector charge and longitudinal component of 
the weak vector current via the Co(V) and Lo(V) multipoles, respectively. The associated 
RMEs, while small, are not negligible-they are about 20 % of the "large" (A) RME in 
3 Si capture, see Table IDG. These 1 S transitions are inhibited by an isospin selection rule, 



indeed they vanish at q=0, since in this limit 



C (q; V) ^ -L ^ *U = ~^=T ± , (6.4) 
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and 



L (q;V) = — 
Q 



H, / dxj (qx)Y o(x) p(x; V) 



1 

<1 



H. 



47T 



(6.5) 



where the expression for L (V) has been obtained by integrating Eq. ( |2.18p by parts, and then 
using the continuity equation to relate V • j(x; V) to the commutator —i[H, p(x; V)]. The 
4 He andp 3 He states have total isospins T, T z =0,0 and 1,1, respectively, ignoring additional, 
but very small, isospin admixtures induced by isospin-symmetry-breaking components of the 
interaction. Therefore matrix elements of the (total) isospin raising or lowering operators 
T± between these T, T z states vanish. 

Equation ( |6.5[) shows that, if the initial and final CHH wave functions were to be exact 
eigenfunctions of the AV18/UIX Hamiltonian, then one would expect, neglecting the kinetic 
energy of the recoiling 4 He: 



E% — Ea 



C (q;V) 



(6.6) 



where E 3 and Ea are the three- and four-nucleon ground-state energies. Note that the 
relation above is in fact valid for any Cj(q; V) and Lj(q;V) multipoles. For g=19.2 MeV/c 
the ratio L /Co is expected to be 1.07, which is in perfect agreement with that obtained in the 
calculation, when the two-body current contributions are taken into account, see Table fTTTi . 
As already discussed in Sec. [IV A| , the present model for the weak vector current satisfies 
current conservation with the Vq part of the nucleon-nucleon interaction (either AV14 or 
AV18). The spin-orbit and quadratic momentum-dependent components of the interaction, 
however, require additional short-range currents that have been neglected in this work. If 
their contributions were to be completely negligible, then the degree of agreement between 
the expected and calculated values for the ratio L /Co would simply reflect the extent to 
which the present variational wave functions are truly exact eigenfunctions of the AV18/UIX 
Hamiltonian. However, the CHH wave function used here gives a ground-state energy of 
—27.9 MeV for 4 He, which is about 400 keV higher than predicted for the AV18/UIX model 
in GFMC calculations 
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In view of these considerations, the perfect agreement referred 
to above may be accidental. 

Finally, the Ci(V) and -^i(V) RMEs interfere destructively in the cross section (see 
discussion at the end of Sec. |II C|) , substantially reducing the 1 Sq channel contribution to 



the S'-factor, see Table XVj. 



B. 3 Si Capture 

The 3 Si capture is induced by the weak axial charge and current, and weak vector current 
operators via the multipoles Ci(A), Li(A), £i(A), and Mi(V). All earlier studies only 
retained the dominant Li(A) and -Eq(A) transitions. However, as is evident from Table [TX|, 
the Mx(V) and especially Ci(A) RMEs are not negligible. Furthermore, the Ci(A) and 
Li(A) RMEs interfere constructively in the cross section, since their signs are opposite. For 
example, neglecting the C*i(A) contribution would produce an S'-factor value 4.94 x 10~ 20 
keV b, 30 % smaller than the 3 Si total result 6.38 x 10" 20 keV b (see Table pCVTD . 
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The destructive interference between the one- and many-body axial current contributions 
in the Iq(A) and -Eq(A) RMEs, first obtained in Refs. [|ll],[2(J, is confirmed in the present 
work. The axial currents associated with A-excitation play a crucial role. The (suppressed) 
one-body contribution comes mostly from transitions involving the D-state components of 
the 3 He and 4 He wave functions, while the many-body contributions are predominantly due 
to transitions connecting the S-state of 3 He to the D-state of 4 He, or viceversa. To clarify 
this point, it is useful to define the one- and two-body densities 

pW(rr) = ( 4 He| - |r, - R jkl \)0^\p s Ee) , (6.7) 



p( 2 \x) 



( 4 He|53a(z-rtf)Oi?b 3 He) 



(6.8) 



where 0\ is the (lowest order) Gamow- Teller operator of Eq. ( 4.12| ) at q=0, and 0\ Z : } is the 



)(2) 



most important A-excitation current associated with diagrams of type (a) in Fig. [3]. These 
densities are normalized such that 



dx p( a \x) = —contribution 



(6.9) 



In Fig. |] we display separately the contributions to p 



due to transitions involving 
the L=0 — >L=0 and L=2 —>L=2 components of the 3 He and 4 He wave functions. Note that 
the L=0 ^L=2 transitions vanish, since the Gamow- Teller operator has no dependence 
on the spatial coordinates in the q=0 limit. The 0—>0 density, while much larger than 
the 2 — > 2 density, consists of positive and negative pieces, which nearly cancel out in the 
integral. Indeed, out of a total integral of 0.19, the — > and 2 — > 2 contributions are, re- 
spectively, 0.02 and 0.17. It is important to reemphasize that in the — > integral the whole 
contribution comes from the mixed symmetry S'-states of the 3 He and 4 He wave functions, 
since the Gamow- Teller operator, in the q = limit, cannot connect their dominant (sym- 



metric) S-states, as already pointed out in Sec. [I A[ This fact has been analitically verified 
using a simplified form for the nuclear wave functions, given by (for 4 He, as an example): 



.crA/ 



U 



a , + u 



trAi 



13. 



i<j 



w 

Kj 



$1 



where $ 4 = det[p 1i,p J, 2 , n t3, n I4] is the spin-isospin Slater determinant, f c,4 {r) 
and u tT ' (r) are central, spin-spin, and tensor correlation functions, respectively. The non- 
central terms in Eq. ( p. 10 ) generate the S'- and D-state components. 

Finally, in Fig. [|, we display both the density functions p^'(x) and p^ 2 \x). The density 
function p( 2 \x), although much smaller than p^'(x), has no zeros, and consequently its 
integral is comparable to that of p^'(x). 

It is interesting to examine the "small" Mi RME induced by the weak vector current. It 
is dominated by the contributions due to two-body currents, which interfere destructively 
with (and, in fact, are much larger in magnitude than) those from one-body currents. This 
matrix element can be approximately related to that occurring in the n 3 He radiative capture 
at thermal neutron energies RTH. Ignoring isospin-symmetry breaking, one has 



(6.10) 



u aA (r) 
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p 3 He) ^^ T+ | T = i ;Mt = ) , (6.11) 



and hence, in a schematic notation, 



( 4 He | e* ■ jt( 7 ) | n 3 He) ~ -±=( 4 He | e* A • jj( 7 ) | T = 1,M T = 0> 



' ( 4 He|el-jt(V)b 3 He) (6.12) 



2C 

where Cq is the Gamow penetration factor, j 2 ( 7 ) is the electromagnetic current, and use has 



T+ , jV(V) 



to jt( 7 ). Note that 



been made of the CVC relation to relate the commutator 

in the first line of Eq. ( |6.12|) the contribution from the T, = 0,0 (l+3)-state has been 
neglected, since the isoscalar magnetic moment of the nucleon is a factor ~ 5 smaller than 
the isovector one, and the dominant two-body electromagnetic currents are isovector. On the 
basis of Eq. (|6.12|) , one would predict n 3 He radiative capture cross sections, at zero energy, 
of 227 fib, 142 fib, and 480 fib with one-body, one- plus two-body, and full currents-the 
latter include the A-excitation currents treated in perturbation theory (PT), which severely 
overestimates their contribution [IT]. The value 480 fib is almost an order of magnitude 
larger than the measured cross section, (55 ± 3) fib [|16| . Ignoring the A contribution, for 
which the PT estimate is known to be unrealistic, the result obtained with one- and two- 
body currents (the model- independent ones of Sec. [IV A| ), 142 fib, is still too large by a 
factor ~ 2.6. However, the approximations made in Eqs. (|6.11| )- (|6.12|) are presumably too 
rough for a reaction as delicate as the n 3 He capture (see discussion in Sec. [TS]). Indeed, 
this process provides a sensitive testing ground for models of interactions and currents. A 
calculation of its cross section with CHH wave functions is currently underway. 

In Table |X| we list the one-body axial current contributions at two values of q, and 
19.2 MeV/c, corresponding to the lowest and highest momentum transfers allowed by the 
p 3 He kinematics. A number of comments are in order. Firstly, the RME associated with 
the Gamow- Teller operator, labelled NR in the table, has a rather strong dependence on q. 
At q=0 this RME is suppressed (see discussion above). When q > 0, however, the next term 
in the expansion of the plane wave in Eq. ( |4.12|) leads to an operator having the structure 



Ti t ± (Ti rf z , which can connect the "large" S- and D-state components of the bound-state wave 
functions. Its contribution, although of order (qR) 2 — 0.02 (R ~ 1.4 fm is the a-particle 
radius), is not negligible. Secondly, the suppression mechanism referred to above also makes 
the relativistic corrections to the Gamow- Teller operator of Eq. ( 4.13 ) relatively important. 



Thirdly, the induced pseudoscalar term, last term in Eq. ( |4.13j ), is purely longitudinal, and 
itself suppressed, since it is proportional to the NR operator. 

In Table [XT| we report the cumulative contributions to the L\(A) and -Ei(A) RMEs at 
q=0 and 19.2 MeV/c. The momentum transfer dependence of the results originates from 
that of the one-body currents, the mesonic and A-excitation current contributions are, in 
fact, very weakly dependent on q. Note that the results obtained by treating the A-currents 
either with the TCO method or in perturbation theory (PT) differ by 1-2 %. This is 
because the NA axial coupling constant g* A is determined by fitting, independently in the 
TCO and PT schemes, the Gamow- Teller matrix element of tritium /5-decay. This procedure 
severely reduces the model dependence of the weak axial current. Finally, we note that, if 
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the quark model value were to be used for g* A (g* A = 6V2/5g A ), the L X (A) (£i(A)) RME at 
g=19.2 MeV/c would have been -0.489 x KT 1 (-0.716 x 10" 1 ) using the TCO method and 
-0.150 x 10" 1 (-0.234 x 10" 1 ) in the PT treatment, respectively. 



C. P-wave Capture 

There are four P-wave capture channels: 3 Po, 1 Pi, 3 Pi, and 3 P2- Note that l Pi and 3 Pi 
are coupled channels (see Sec. Ill C ). The 3 Pq capture is induced by the weak axial charge 



and the longitudinal component of the weak axial current via the Co (A) and Lo(A) mufti 
poles, respectively. The associated RMEs, as defined in Eq. (|6.1|), are listed in Table |X11 



The two-body axial charge operators of Sec. [IV D| , among which the pion-exchange term is 



dominant, give a ~ 20 % correction to the one-body contribution in the Co (A) RME. The 
Lo (A) RME is about 40 % of, and has the same sign as, the Co (A) RME. This positive 
relative sign produces a destructive interference between these RMEs in the cross section, 
substantially reducing the 3 Pq overall contribution to the 5-factor, as discussed in Sec. [11 C 



The C (A) and L (A) RMEs are expected to have the same sign, as justified by the following 
argument. The C (q; A) multipole operator can be written, in the q — > limit, as 

1 9a 



C (q; A) -> 7=77^- Y] K± °V > P 



i 9a 
4T 2 



^ [ri,± cr i ■ r j , H , (6.13) 

i 

where we have used the approximate relation pj ~ — i m [r» , H ] (violated by the momentum- 
dependent components of the two-nucleon interaction), and in the second line of Eq. (|6.13| ) 
have ignored, in a rather cavalier fashion, terms like n,± en ■ Hri — YiH ■ (TiT% t ±. For the 
L (q; A) multipole we find in the same limit 

i 

and therefore we would expect the Co (A) and Lq(A) RMEs to be approximately in the ratio 



L (A) 2 q 

which, given the rather severe approximations made in deriving Eq. ( |6.13| ), is reasonably 
close to the (one-body) value obtained in the calculation (1.6 versus 2.0). 

The x Pi and 3 Pi captures are induced by the weak vector charge and current, and 
weak axial current via the multipoles Ci(V), ^i(V), Ex(V), and Mi (A). The calculated 
values for the associated RMEs are listed in Tables [XIII| and [XIV[ The RME magnitudes 



of the weak vector transitions in 3 Pi capture are much smaller than those in : Pi capture. 
In the long- wavelength approximation, the one-body Ci(V), L X (V), and -Ei(V) multipoles 
are independent of spin, and therefore cannot connect the dominant part of the 3 Pi wave 
function, which has total spin S—l, to the S-wave component of 4 He, which has S=0. This is 
not the case for the x Pi channel, in which the total spin S=0 term is in fact largest. Indeed, 
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because of this suppression, the two-body weak vector charge and current contributions are 
found to be dominant in 3 Px capture. The situation is reversed for the axial transition, since 
there the spin-flip nature of the Mi (A) multipole makes the associated RME in 3 P! larger 
than that in 1 P 1 (in absolute value). 

The Ei(V) operator can be shown to have the long- wavelength form 



E 1 (q;V) = -^[H,C 1 (q;V)} , (6.16) 
Q 

and so the -Ei(V) and Ci(V) RMEs would be expected to be in the ratio 

— — — -~V2 -1.51, 6.17 

Cx(q;V) q 

assuming the validity of the long-wavelength approximation, and that the CHH wave func- 
tions are truly exact eigenfunctions of the Hamiltonian. We reiterate here that the currents 
used in the present work satisfy the continuity equation only with the vq part of the AV14 

and AV18 interactions, namely in momentum space q • j(q; V) = T + t> 6 , pj^q; V) . The 

currents induced by the momentum dependent components of the interactions, such as the 
spin-orbit term, have been neglected. Thus the ratio obtained in the calculation is 1.34 
for the x Pi channel, somewhat smaller than the expected value presumably because of the 
"missing" currents and the approximate eigenstate property satisfied by the present CHH 
(variational) wave functions. These same cautionary remarks also apply to the comparison 
between the Ci(V) and £i(V) RMEs, which should be related to each other via Eq. (|6.6| ). 

The situation is more delicate in 3 Pi capture, since this transition is suppressed. Here 
the long-wavelength approximation of the -E'i(V) multipole is inadequate, and higher order 
terms in the power expansion in q need to be retained, so called retardation terms. In fact 
the situation is closely related to that of electric dipole transitions in pd radiative capture 



at very low energies (0-100 keV). We refer the reader to Ref. |T8| for a thorough discussion 
of these issues. 

The 3 P2 capture is induced by the weak axial charge and current, and weak vector current 
operators via the multipoles C 2 (A), L 2 (A), E 2 (A), and M 2 (V). The associated RMEs are 
listed in Table |XV[ The L 2 (A) and E 2 (A) RMEs are comparable to the Li(A) and E^A) 



RMEs in 3 Si capture, and are dominated by the contributions of one-body currents. In 
fact, the latter can now connect the large S-wave components of the three- and four-nucleon 
bound states. The density function p^^x), defined in analogy to Eq. ( |6 . T[ ) (but for the 3 P 2 
channel), is displayed in Fig. [| and should be compared to that in Fig. [| for 3 Si capture. 
While smaller in magnitude than the latter-after all, the 3 P 2 transition is inhibited with 
respect to the 3 Si transition by a factor ~ qR and the presence of the centrifugal barrier-the 
3 P 2 density has the same sign, and therefore its integral turns out to be comparable to that 
of the 3 Si density. 



D. Model Dependence 

In Table [XVI] we list, for all S- and P-wave channels, the S-factor values obtained with 
the AV18/UIX, AV18, and AV14/UVIII interactions. Note that the sum of the channel 
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contributions is a few % smaller than the total result reported at the bottom of the table 
(see end of Sec. [II C|) . The iVA axial coupling constant is determined by fitting the Gamow- 
Teller matrix element in tritium /3-decay, within each given Hamiltonian model. As a result 
of this procedure the model dependence of the S-factor predictions is substantially reduced. 

Inspection of Table |XV1| shows that inclusion of the three-nucleon interaction reduces 
the total S'-factor by about 20 % (compare the AV18 and AV18/UIX results). This decrease 
is mostly in the 3 Si contribution, and can be traced back to a corresponding reduction in 
the magnitude of the one-body axial current matrix elements. The latter are sensitive to the 
triplet scattering length, for which the AV18 and AV18/UIX models predict, respectively, 
10.0 fm and 9.13 fm (see Table |TTT|) . 

The comparison between the AV18/UIX and AV14/UVIII models, which both reproduce 
the measured bound-state properties and low-energy scattering parameters of the three- and 
four-nucleon systems, suggests a rather weak model dependence. It is important to reiterate 
that this is accomplished by virtue of the procedure used to constrain the axial current. 
Indeed, the AV18/UIX and AV14/UVIII 3 Si contributions to the S'-factor obtained with 
one-body currents only are, respectively, 26.4 x 10 -20 keV b and 35.8 x 10~ 20 keV b. This 
difference is presumably due to the stronger tensor component of AV14 as compared to that 
of AV18. 

Finally, the 3 Si contribution to the S-factor obtained with the AV14/UVIII model in the 
present work, 6.60 x 10~ 20 keV b, is to be compared with the older prediction of Ref. [11 1, 
1.3 x 10~ 20 keV b. It is important to point out that the older calculation (i) used the 
long- wavelength form of the E\(A) and L\(A) operators, (ii) ignored the contributions of 
transitions induced by the axial charge and vector current, (iii) retained only the leading 
non-relativistic (Gamow- Teller) term of the single-nucleon axial current, and (iv) employed 
bound and continuum wave functions, obtained with the Variational Monte Carlo (VMC) 
method. In regard to this last point, we note that, for example, the Ei(q = 0; A) RME 
calculated in Ref. [O with the Gamow- Teller operator is 0.613 x lO^ 1 fm 3//2 versus a value 
of 0.119 fm 3 / 2 obtained here. The factor ~ 2 increase is only due to differences in the wave 
functions. The present CHH wave functions are expected to be more accurate than the 



VMC wave functions of Ref. 1 1 
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TABLES 



TABLE I. The hep S-factor, in units of 10 keV b, calculated with CHH wave functions 
corresponding to the AV18/UIX Hamiltonian model, at j> 3 He cm. energies E=0, 5, and 10 keV. 
The rows labelled "one-body" and "full" list the contributions obtained by retaining the one-body 
only and both one- and many-body terms in the nuclear weak current. The contributions due the 
3 Si channel only and all S- and P-wave channels are listed separately. 





E= 


=0 keV 


E= 


=5 keV 


E= 


10 keV 






S+P 




S+P 


3 Si 


S+P 


one-body 


26.4 


29.0 


25.9 


28.7 


26.2 


29.3 


full 


6.38 


9.64 


6.20 


9.70 


6.36 


10.1 



TABLE II. Binding energies in MeV of 4 He calculated with the CHH method using the AV18 
and AV18/UIX, and the older AV14 and AV14/UVIII, Hamiltonian models. Also listed are the 
corresponding "exact" GFMC results p2| , |33[ and the experimental value. 



Model 




CHH 




GFMC 


AV18 

AV18/UIX 
AV14 

AV14/UVIII 




24.01 
27.89 
23.98 
27.50 




24.1(1) 
28.3(1) 
24.2(2) 
28.3(2) 


EXP 








28.3 




TABLE III. Binding energies, B3, of 3 He, and p 3 He singlet and triplet S-wave scattering 
lengths, a s and at, calculated with the CHH method using using the AV18 and AV18/UIX, and 
the older AV14 and AV14/UVIII, Hamiltonian models. The corresponding experimental values are 
also listed. 


Model 


S 3 (MeV) 


a s (fm) 




a t (fm) 




AV14 
AV18 

AV14/UVIII 
AV18/UIX 


7.03 
6.93 
7.73 
7.74 


12.9 
11.5 




10.0 
9.24 
9.13 




EXP 


7.72 


10.8+2.6 | 


41 


] 8.1+0.5 [ 
10.2+1.5 


1 

271 
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TABLE IV. Phase-shift and mixing-angle parameters (in deg) for p 3 He elastic scattering at 
cm. energy of 1.2 MeV, calculated with the CHH method using the AV18 and AV18/UIX Hamil- 
tonian models. The corresponding experimental values obtained in the phase-shift analysis of 



Ref. [41 1 are also listed. 



parameter AV18 AV18/UIX PSA 

% -33.3 -31.3 -27.4±3.5 

3 Si -28.8 -27.1 -26.5±0.6 

3 P n 4.1 3.2 2.6±0.6 



o 

6.5 5.5 4.2±1.5 

e(l-) -14.7 -13.2 -7.8±0.6 



3 Pi 8.1 7.4 10.1±0.5 

3 P 2 7.7 6.9 8.9±0.5 



TABLE V. Contributions to the Gamow- Teller (GT) matrix element of tritium /3-decay, ob- 
tained with the CHH trinucleon wave functions corresponding to the AV18/UIX Hamiltonian 
model. The rows labelled "one-body NR"and "one-body RC'list the contributions associated with 
the single- nucleon axial current operators of Eq. (4.12) and Eq. ( |4.13| ), respectively, while the row 



labelled "mesonic" lists the sum of the contributions due to the ir-, p-, and p7r-exchange axial cur- 
rent operators of Eqs. fl4.32|) -( fl.34| ) with cutoff masses 

K = 4.8 fur 1 and A p = 6.8 fen" 1 . Finally, 
the rows labelled "A-g^", u A-g~ A ", and "A-renormalization"list, respectively, the contributions 
associated with panels (a)-(b), (c)-(d) and (f), and (e) and (g)-(j), of Fig. |3|. The cumulative 
result reproduces the "experimental value" 0.957 for the GT matrix element |^|, once the change 
in normalization of the wave functions due to the presence of A-components is taken into account. 

GT matrix element 

one-body NR 0.9218 

one-body RC -0.0084 

mesonic 0.0050 

A-g* A 0.0509 

A-g A 0.0028 

A-renormalization 0.0074 



TABLE VI. The values of the N — > A axial coupling constant g* A in units of gA, when the 
A-isobar degrees of freedom are treated in perturbation theory (PT), or in the context of a TCO 
calculation based on the AV28Q interaction. The purely nucleonic CHH wave functions correspond 
to the AV18/UIX Hamiltonian model. 



A-isobar treatment 


9* A /9A 


PT 


1.224 


TCO 


2.868 
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TABLE VII. The wave function normalization ratios (^n+a | ^ n+a) / | obtained for the 
bound three- and four-nucleon systems, when the TCO calculation is based on the AV28Q inter- 
action. The purely nucleonic CHH wave functions | ^) correspond to the AV18/UIX Hamiltonian 
model. 



Model 


3 H 


3 He 


4 He 


AV28Q 


1.0238 


1.0234 


1.0650 



TABLE VIII. Cumulative contributions to the reduced matrix elements (RMEs) Co(q; V) and 
Lo(q; V) in 1 So capture at zero p 3 He cm. energy. The momentum transfer q is 19.2 MeV/c, and 
the results correspond to the AV18/UIX Hamiltonian model. The row labelled "one-body" lists 
the contributions associated with the operators in Eq. (4.5) for the weak vector charge p(V) and 
Eq. (4.8) for the weak vector current j(V); the row labelled "mesonic" lists the results obtained 
by including, in addition, the contributions associated with the operators in Eqs. (4.30)-(4.31) for 
p(V), and Eqs. (4.16)-(4.17) for j(V). The A terms in p(V) are neglected, while those in j(V) are 
purely transverse and therefore do not contribute to the Lq RME. Note that the RMEs are purely 
real and in fm 3 / 2 units. 











C (g;V) 


L (q; V) 


one-body 


-0.857 x 10~ 2 


-0.864 x 10" 2 


mesonic 


-0.856 x 10" 2 


-0.919 x 10~ 2 



TABLE IX. Cumulative contributions to the reduced matrix elements (RMEs) Ci(g;A), 
Li(q; A), E\{q; A) and Mi(q; V) in 3 Si capture at zero p 3 He cm. energy. The momentum trans- 
fer q is 19.2 MeV/c, and the results correspond to the AV18/UIX Hamiltonian model. The row 
labelled "one-body" lists the contributions associated with the operators in Eq. (4.10) for the weak 
axial charge p(A), Eq. (4.11) for the weak axial current j(A), and Eq. (4.8) for the weak vector 
current j(V); the row labelled "mesonic" lists the results obtained by including, in addition, the 
contributions associated with the operators in Eqs. (4.35)-(4.37) for p(A), Eqs. (4.32)-(4.34) for 
j(A), and Eqs. (4.16)-(4.17) for j(V); finally, the row labelled "A"lists the results obtained by also 
including the contributions of the operators in Eqs. (4.50)-(4.51) for p(A), Eqs. (4.48)-(4.49) for 
j(A), and Eqs. (4.52)-(4.53) for j(V). The A contributions in both p(A) and j(A) are calculated 
with the TCO method, and take into account the change in normalization of the wave functions 
due to the presence of A-components. Those in j(V) are calculated in perturbation theory. Note 
that the RMEs are purely imaginary and in fm 3 / 2 units. 





Cife A) 


Li(g;A) 


E 1 (q;A) 


Mi(g;V) 


one-body 


0.147 x 10" 1 


-0.730 x 10" 1 


-0.106 


0.333 x 10" 2 


mesonic 


0.156 x 10" 1 


-0.679 x lO" 1 


-0.984 x 10" 1 


-0.263 x 10~ 2 


A 


0.155 x 10- 1 


-0.293 x 10- 1 


-0.440 x 10- 1 


-0.484 x 10~ 2 
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TABLE X. One-body contributions, at momentum transfers q=0 and 19.2 MeV/c, to the re- 
duced matrix elements (RMEs) L±(q; A) and Ei(q; A) in 3 Si capture at zero p 3 He cm. energy. The 
results correspond to the AV18/UIX Hamiltonian model. The rows labelled "NR"and "RC'list the 
contributions obtained with the operators of Eq. (4.12) and Eq. (4.13), respectively; the row la- 
belled "IPS"lists the contribution of the induced pseudoscalar current only (last term in Eq. (4.13)). 
Note that the RMEs are purely imaginary and in fm 3 / 2 units. 





Li(g; A) 




Ei(q;A) 






5=0 MeV/c 


5=19.2 MeV/c 


5=0 MeV/c 


5=19.2 MeV/c 


NR 


-0.726 x KT 1 


-0.586 x lO" 1 


-0.103 


-0.838 x lO" 1 


RC 


-0.154 x 1Q- 1 


-0.145 x lO" 1 


-0.220 x 1Q- 1 


-0.219 x 1Q- 1 


IPS 




0.741 x 10~ 3 







TABLE XI. Cumulative contributions, at momentum transfers q=0 and 19.2 MeV/c, to the 
reduced matrix elements (RMEs) L\(q;A) and E\{q\A) of the weak axial currrent in 3 Si capture 
at zero p 3 He cm. energy. The results correspond to the AV18/UIX Hamiltonian model. The 
row labelled "one-body" lists the contributions associated with the operator in Eq. (4.11); the row 
labelled "mesonic" lists the results obtained by including, in addition, the contributions associated 
with the operators in Eqs. (4.32)-(4.34); finally, the rows labelled "A-TCO"and A-PT list the 
results obtained by also including the contributions of the operators in Eqs. (4.48)-(4.49), calculated 
either in the TCO scheme or in perturbation theory (PT). The A-TCO results also take into account 
the change in normalization of the wave functions due to the presence of A-components. Note that 



the RMEs 


are purely imaginary and in 


fm 3 / 2 units. 








17(5; A) 






Ei(q;A) 




5=0 MeV/c < 


? =19.2 MeV/c 


5=0 MeV/c 


5=19.2 MeV/c 


one-body 


-0.880 x lO" 1 


-0.730 x 10" 1 


-0.125 


-0.106 


mesonic 


-0.829 x 10" 1 


-0.679 x 10" 1 


-0.117 


-0.984 x 10- 1 


A-TCO 


-0.440 x 10- 1 


-0.293 x 10" 1 


-0.625 x 10" 


1 -0.440 x 10- 1 


A-PT 


-0.447 x 1Q- 1 


-0.298 x 1Q- 1 


-0.631 x 10- 


1 -0.443 x 10" 1 
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TABLE XII. Cumulative contributions to the reduced matrix elements (RMEs) Co(q; A) and 
Lo(q; A) in 3 Po capture at zero p 3 He cm. energy. The momentum transfer q is 19.2 MeV/c, and 
the results correspond to the AV18/UIX Hamiltonian model. The row labelled "one-body" lists 
the contributions associated with the operators in Eq. (4.10) for the weak axial charge p(A) and 
Eq. (4.11) for the weak axial current j(A); the row labelled "mesonic" lists the results obtained 
by including, in addition, the contributions associated with the operators in Eqs. (4.35)-(4.37) for 
p(A), and Eqs. (4.32)— (4.34) for j(A); finally, the row labelled "A"lists the results obtained by also 
including the contributions of the operators in Eqs. (4.50)— (4.51) for p(A), and Eqs. (4.48)— (4.49) 
for j(A). The A contributions in both p(A) and j(A) are calculated with the TCO method, 
and take into account the change in normalization of the wave functions due to the presence of 
A-components. Note that the RMEs are purely imaginary and in fm 3 / 2 units. 





Cb(</; A) 




one-body 


0.371 x 10- 1 


0.182 x 10" 1 


mesonic 


0.444 x lO" 1 


0.183 x 10- 1 


A 


0.459 x lO" 1 


0.188 x 1Q- 1 



TABLE XIII. Cumulative contributions to the reduced matrix elements (RMEs) Ci(g;V), 
Li(q;V), Ei(q;Y) and Mi(g;A) in 1 Pi capture at zerop 3 He cm. energy. The momentum trans- 
fer q is 19.2 MeV/c, and the results correspond to the AV18/UIX Hamiltonian model. The row 
labelled "one-body" lists the contributions associated with the operators in Eq. (4.5) for the weak 
vector charge p(V), Eq. (4.8) for the weak vector current j(V), and Eq. (4.11) for the weak axial 
current j(A); the row labelled "mesonic" lists the results obtained by including, in addition, the 
contributions associated with the operators in Eqs. (4.30)-(4.31) for p(V), Eqs. (4.16)-(4.17) for 
j(V), and Eqs. (4.32)-(4.34) for j(A); finally, the row labelled "A" lists the results obtained by also 
including the contributions of the operators in Eqs. (4.52)-(4.53) for j(V), and Eqs. (4.48)-(4.49) 
for j(A). The A contributions in j(A) are calculated with the TCO method, and take into account 
the change in normalization of the wave functions due to the presence of A-components. Those 
in j(V) are calculated in perturbation theory. The A terms in p(V) are neglected, while those in 
j(V) are purely transverse and therefore do not contribute to the L\ RME. Note that the RMEs 
are purely real and in fm 3 / 2 units. 





Ci(g;V) 


Li(q;V) 


Ei(q;V) 


M\(q; A) 


one-body 


-0.222 x 10- 1 


-0.162 x 10- 1 


-0.231 x 10- 1 


-0.100 x 10" 2 


mesonic 


-0.222 x 10- 1 


-0.209 x 10- 1 


-0.298 x 10- 1 


-0.779 x 10~ 3 


A 






-0.298 x 1Q- 1 


-0.809 x 10~ 3 
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TABLE XIV. Cumulative contributions to the reduced matrix elements (RMEs) Ci(q;V), 
E\(q;\) and M\(q;A) in the 3 Pi capture at zero p 3 He cm. energy. The momentum 
transfer q is 19.2 MeV/c, and the results correspond to the AV18/UIX Hamiltonian model. Nota- 
tion as in Table KIH| . Note that the RMEs are purely real and in fm 3 / 2 units. 

CT(g;V) gfaivj ^(g;V) MT(g;A) 

one-body 0.953 x 10~ 3 0.118 x 10~ 2 0.521 x 10~ 3 0.304 x 10" 1 

mesonic 0.217 x 10~ 2 0.174 x 10~ 2 0.128 x 10~ 2 0.304 x 10" 1 

A 0.127 x 10~ 2 0.303 x 10" 1 



TABLE XV. Cumulative contributions to the reduced matrix elements (RMEs) C2(q;A), 
L%(q;A), E2(q;A) and M2(<?;V) in the 3 P2 capture at zero p 3 He cm. energy. The momentum 
transfer q is 19.2 MeV/c, and the results correspond to the AV18/UIX Hamiltonian model. Nota- 
tion as in Table [[XL Note that the RMEs are purely imaginary and in fm 3//2 units. 

QfajA) T 2 (q;A) Ej(q;A) M(q;y) 

one-body -0.146 x 10~ 3 0.236 x 10" 1 0.292 x 10" 1 -0.110 x 10" 2 

mesonic -0.114 x 10~ 3 0.236 x 10" 1 0.293 x 10" 1 -0.116 x 10~ 2 

A -0.988 x 10" 4 0.238 x 10" 1 0.295 x 10" 1 -0.118 x 10~ 2 



TABLE XVI. Contributions of the S- and P-wave capture channels to the hep S'-factor at zero 
p 3 He cm. energy in 10~ 20 keV b. The results correspond to the AV18/UIX, AV18 and AV14/UVIII 



Hamiltonian models. 





AV18/UIX 


AV18 


AV14/UVIII 




0.02 


0.01 


0.01 


3 S! 


6.38 


7.69 


6.60 


3 Po 


0.82 


0.89 


0.79 


X Pl 


1.00 


1.14 


1.05 


3 Pl 


0.30 


0.52 


0.38 


3 P 2 


0.97 


1.78 


1.24 


TOTAL 


9.64 


12.1 


10.1 
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FIGURES 
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FIG. 1. Electron energy spectrum for the ratio between the Super-Kamiokande 825-days data 



and the expectation based on unoscillated B neutrinos |lCj]. The data were extracted graphically 
from Fig. 8 of Ref. Q . The 5 curves correspond respectively to no hep contribution (dotted line) , 
and an enhancement a of 2.2 (solid red line), 4.4 (solid blue line), 10 (solid green line) and 20 
(long-dashed green line). 
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FIG. 2. Differential cross section cr(0) as function of the cm. scattering angle 9 at cm. energy 
of 1.2 MeV. The experimental data are taken from Ref. |j62"[ . The long-dashed and solid lines 
correspond, respectively, to the CHH calculations with the AV18 and AV18/UIX Hamiltonian 
models. 
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FIG. 3. Diagrammatic representation of the operators included in O(A), due to the one-body 
current and charge operators 0^(N — » A), 0^(A — > N), and 0^(A — » A), given in 
Eqs. (|4.48| )- (|4.53| ), and to the transition correlations U AN , U N , U^ A , and corresponding her- 

mitian conjugates. Thin, thick, and dashed lines denote, respectively, nucleons, A isobars, and 

' 't , 

transition correlations U BB or U BB , with B,B = N, A. 
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FIG. 4. Contributions to the density function p^\x), defined in Eq. fl6.7|), due to transitions 
involving the L = — ► L = (filled circles) and L = 2 — > L = 2 (opaque squares) components in the 
3 He and 4 He wave functions. Note that the 2 — ► 2 density function has been multiplied by a factor 
of 10, for ease of presentation. 
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FIG. 5. Density functions (filled circles) and p( 2 )(a;) (opaque squares), denned in 

Eqs. 0-(PD- 
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FIG. 6. Density function p^(x) defined in Eq. (|6.7|) in the 3 P2 capture channel. 



54 



